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The radiations accompanying the decay of Mn® and of 
Co* were studied by means of the magnetic lens spec- 
trometer and coincidence methods. Both of these nuclei 
decay to stable Fe® as the product nucleus. Four excited 
states of the Fe® nucleus are definitely established with 
excitation energies 0.845 Mev, 2.11 Mev, 2.66 Mev, and 
2.98 Mev, respectively. De-excitation of the three high 
energy states always leads to the 0.845-Mev level and 
thence to the ground state. The negatron spectrum of 
Mn*® consists of three groups of maximum energies 2.86 
Mev, 1.05 Mev, and 0.73 Mev, respectively, and relative 
abundance 60:25:15 leading to the first, third, and 
fourth of the excited states mentioned above. The positron 
spectrum of Co consists mainly of a single group of 
maximum energy 1.50 Mev, leaving the Fe* nucleus in 
the 2.11-Mev state. Orbital electron capture also takes 
place involving several other excited states. Altogether 


the energies of eight gamma-rays emitted by the Fe™* 
nucleus are determined: 0.845 Mev, 1.26 Mev, 1.74 Mev, 
1.81 Mev, 2.01 Mev, 2.13 Mev, 2.55 Mev, 3.25 Mev. 
Within the experimental uncertainty five of these energies 
appear to be the multiples 2, 3, 4, 5, and 6 of 0.425 Mev. 
All of the beta-ray spectra involved seem to have the 
“allowed” shape of the Fermi theory. The significance of 
this and of the relative probability of orbital electron 
capture is discussed in terms of the tensor interaction. 
It appears that all of the transitions should involve nuclear 
angular momentum changes AJ= +1 or AJ=0, with or 
without parity change. From the disintegration schemes 
one obtains the mass differences between neutral atoms, 
namely, Mn*—Fe*=3.98 mMU and Co*—Fe*=4.96 
mMU, and the threshold for the reaction Fe*(p, n)Co™: 
5.47 Mev. 





INTRODUCTION 


HE preceding papers of this series' have 
described investigations leading to the 
assignment of energy levels in product nuclei of 
radioactive disintegrations. The question natu- 
rally arises whether any particular radioactive 
decay involves all, or even most, of the energeti- 
cally available levels of the product nucleus. 
The results of the previous investigations indi- 
cate that only a few of the possible transitions 
take: place with sufficient frequency to be de- 
tected by our present methods. Thus, for ex- 
*Now with the Radiological Committee, Canadian 
National Research Council, Montreal, Canada. 
'Phys. Rev. 60, 470L (1941); 60, 544 (1941); 61, 


686 (1942); 62, 3 (1942) ; 64, 268 (1943), hereafter referred 
to as papers I to V. 


ample in the decay of I'*° (paper V) only two of 
five possible beta-ray transitions and four of ten 
possible gamma-ray transitions between the 
known levels are actually observed. In the case 
of Br® (paper II) only one of four possible 
beta-ray transitions and three of six gamma-ray 
transitions are found. So far it has always been 
possible to account for the absence of the other 
transitions on the basis of known selection rules 
for gamma-ray transitions and the Gamow- 
Teller selection rules for beta-ray emission. 
Consistent assignments of total angular momenta 
and parity have been possible, in this way, 
although these assignments are by no means 
unique. Generally the average level spacings are 
found to be of the same order of magnitude as 
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Fic. 1. The beta-ray spectrum of Mn**. The open 
circles and heavy line represent the experimental data. 
The full circles indicate the three partial spectra, as derived 
from the Fermi plot. 


those found in transmutation? and _ inelastic 
scattering’ experiments. A comparison of the 
results of heavy particle studies with radioactive 
disintegration schemes should permit an estimate 
of the fraction of all states likely to be detected 
by any one of the methods. Unfortunately very 
few nuclei have been studied with sufficient 
accuracy by more than one of these methods 
giving information about low lying excited states. 

It will probably be necessary to improve the 
resolution of the heavy particle studies, possibly 
by substituting magnetic analysis for the present 
absorption methods before they can be used for 
actual identification of energy levels. It is doubt- 
ful whether, e.g., the accuracy of any of the 
present transmutation and scattering experi- 
ments would be sufficient to resolve the two 
known excited states of Co*® (paper IV) or to 
decide the order of the two states in Ni®.‘ 

One may hope, however, to obtain valuable 
information by studying the radioactive decay 
of different parent nuclei to the same product 
nucleus. This may occur either by the decay of 
independent (genetically unrelated) isomers or 
by the decay of radioactive isobars differing by 
two in atomic number. The first case occurs 
rather rarely such as in Co® which is, however, 
difficult to study**® because of the short life of 
one of the isomeric activities. The convergent 
decay of two isobars occurs fairly frequently, 

? E. F. Shrader and E. Pollard, Phys. Rev. 59, 277 (1941). 


3R. H. Dicke and J. Marshall, Jr., Phys. Rev. 63, 86 
(1943). 

4M. Deutsch and L. G. Elliott, Phys. Rev. 62, 558A 
(1942). 

5H. E. Nelson, M. L. Pool, and J. D. Kurbatov, Phys. 
Rev. 62, 1 (1942). 


but only a few of these isobaric pairs of activities 
are suitable for our purpose. In some cases (e.g., 
Cr®) the half-life of one or both of the activities 
is inconveniently short. Others (e.g., Cu®) lack 
interest because one or both of the disintegrations 
involves only the ground state of the product 
nucleus. 

Preliminary investigations indicated that the 
isobaric pair Mn**—Co* is very well suited for 
a study of the levels of Fe®® since both disinte- 
grations involve several gamma-rays and proceed 
sufficiently slowly for careful investigation. A 
brief report on some results of such a study was 
made in a Letter to the Editor® and the present 
paper extends these results and describes the 
experiments leading to them. 


A. Mn** 


The 2.59-hr. activity of manganese is definitely 
assigned to the species Mn* since it is produced 
by slow neutron bombardment of stable Mn‘, 
The radiations accompanying its decay have 
been studied by a number of authors. We shall 
refer here only to some of the more recent results, 
Townsend’ found that the beta-ray spectrum, 
observed in a 180° focusing spectrometer, appears 
to consist of two groups with maximum energies 
2.88 Mev and 1.01 Mev, respectively. This is in 
good agreement with the results of Langer, 
Mitchell, and McDaniel* and of Dunworth® who 
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Fic. 2. Fermi plot of the spectrum shown in Fig. 1. 
The end points of the low energy spectra deduced from 
this plot are in excellent agreement with the level scheme 
found from gamma-ray energy and coincidence measure- 
ments. 


‘L. G. Elliott and M. Deutsch, Phys. Rev. 63, 321L 
(1943). 

7A. A. Townsend, Roy. Soc. Proc. A177, 362 (1941). 

§L. M. Langer, A. C. G. Mitchell, and P. W. McDaniel, 
Phys. Rev. 56, 422 (1939). 

* J. V. Dunworth, Nature 143, 1065 (1939). 
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found that all beta-rays of energy greater than 
about 1 Mev are accompanied, on the average, 
by the same gamma-rays whose energy the latter 
author estimated to be about 0.60 Mev, while 
the beta-rays of lower energy are accompanied 
by harder gamma-rays. Langer et al.* estimated 
the energy of these harder gamma-rays as 1.7 
Mev. Gamma-gamma coincidences showed that 
some gamma-rays are emitted in cascade, and 
Dunworth® proposed a disintegration scheme in 
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Fic. 3. Secondary electron spectra due to the gamma- 
rays of Mn**. The two curves do not correspond to sources 
of equal strength. 


which 60 percent of the beta-rays are accom- 
panied by a single 0.6-Mev gamma-ray, while the 
remainder is accompanied by a 1.7-Mev and a 
0.6-Mev gamma-ray in cascade. However, Cur- 
ran, Dee, and Strothers'® found two gamma-rays 
of 0.91-Mev and 2.01-Mev energy, respectively, 
by studying the Compton electrons ejected by 
the gamma-rays of Mn*® in a 180° focusing 
spectrometer. Deutsch and Roberts" studying 
photoelectrons as well as Compton recoils found 
values of 0.835 and 2.13 Mev, respectively. 
Although these measurements were not incom- 
patible with Dunworth’s disintegration scheme, 
the disagreement between the best beta-ray’ and 
gamma-ray" energy determinations was serious 
enough to make the proposed scheme rather 
unsatisfactory. Recent improvements in spec- 
trometer technique have allowed us to discover 
a third gamma-ray of energy 1.81 Mev which 
resolves these difficulties and leads to a satis- 
factory modification of the disintegration scheme. 


1S. C. Curran, P. I. Dee, and J. E. Strothers, Proc. 
Roy. Soc. A174, 546 (1940). 

"M. Deutsch and A. Roberts, Phys. Rev. 60, 362L 
(1941). 


PREPARATION OF SOURCES 


Sources of Mn*® were usually prepared by 
exposing a saturated solution of KMnQ, to slow 
neutrons from the M.I.T. cyclotron for several 
hours. The solution was then filtered through a 
fritted glass filter, and most of the activity 
was found in the MnO, precipitate. Filtering 
the KMnQ, solution through paper increases the 
amount of precipitate without increasing the 
activity recovered, probably because of the re- 
ducing action of the paper. 

Thin sources for beta-ray studies were pre- 
pared by depositing a drop of active MnCl, 
solution and then a drop of (NH,).S solution on 
a thin flake of mica and evaporating to dryness. 
The source was then heated for several minutes 
to drive off excess reagent. This procedure is 
preferable to simple evaporation of a chloride 
solution since the finely divided sulfide precipi- 
tate forms a much more uniform and coherent 
deposit than a salt which crystallizes from the 
solution. The mica was mounted in the spec- 
trometer without backing to minimize back 
scattering of the beta-rays. 

For the study of gamma-ray energies, active 
MnO, was sealed into a brass capsule of sufficient 
wall thickness to stop all primary beta-rays. 
Suitable lead and bismuth radiators were fixed 
to this capsule with wax for the measurement of 
photoelectron distributions. 


THE BETA-RAY SPECTRUM OF Mn* 


Figure 1 shows the momentum distribution of 
the beta-rays of Mn**, obtained in the- magnetic 
lens spectrometer. We believe that source and 
backing were thin enough to meet the rigorous 
criteria for reliable beta-ray measurements. 
Figure 2 shows the same data as a Fermi plot. 
Our results agree very closely with those obtained 
by Townsend.’ It is quite obvious from Fig. 2 
that there are at least two groups of beta-rays 
with maximum energies 2.86 Mev and about 
1 Mev, respectively. The particular analysis of 
the Fermi plot into three straight lines as shown 
in Fig. 2 is based on evidence about the disinte- 
gration scheme obtained from the gamma-ray 
energies and coincidence measurements. It is by 
no means to be taken for granted that the Fermi 
plot can be analyzed into straight lines repre- 
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senting various groups since some simple ‘‘for- 
bidden”’ spectra are known not to be represented 
by straight Fermi plots,” and it is unlikely that 
all of the groups emitted by Mn* are emitted in 
“‘allowed”’ transitions. The significance of the 
fact that, nevertheless, when the beta-ray spec- 
trum of Mn*® is analyzed into groups with the 
maximum energies predicted by independent 
information about the disintegration scheme, 
each group is represented by a straight Fermi 
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Fic. 4. a. Disintegration scheme of Mn**. b. Level 
scheme for the most abundant mode of disintegration of 
Co**. c. Proposed level scheme for the transitions to Fe**. 
Solid lines and arrows indicate levels and transitions 
which appear well established. Dotted lines and arrows 
indicate an arrangement of the observed gamma-rays of 
Co**, consistent with the probable selection rules. 











plot, will be discussed further in Section C 
below. This result is similar to those obtained in 
the case of Fe®® (paper IV) and of I'*° (paper V). 
The three groups are shown in Fig. 1. Their 
maximum energies are 0.75+0.1 Mev, 1.05+0.03 
Mev, and 2.86+0.05 Mev, respectively, and 
their relative abundance 15 : 25 : 60. 


THE GAMMA-RAY SPECTRUM 


Figure 3 shows the spectra of secondary 
electrons produced in lead radiators of two 
different thicknesses (and in the brass capsule 


2E. J. Konopinski and G. E. Uhlenbeck, Phys. Rev. 
60, 308 (1941). 


containing the source) by the gamma-rays of 
Mn**, The observed counting rates were not 
divided by Hp as is customary, since the heights 
of the photoelectron lines, representing intensities 
of nearly monokinetic groups of electrons, should 
not be so divided."* There appear three pro- 
nounced photoelectron lines due to K electrons 
of lead ejected by gamma-rays of 0.845+0.015- 
Mev, 1.81+0.04-Mev, and 2.13+0.04-Mev 
energy, respectively. The incompletely resolved 
L-electron group due to the lowest energy 
gamma-ray gives rise to the peak on the high 
energy side of the corresponding K electrons. 
The continuous background consists of Compton 
recoil electrons due to the same three gamma- 
rays. Only the energy range above 0.7 Mev is 
shown since no gamma-rays were found below 
that energy. 

The resolution obtained in the experiments of 
Curran, Dee, and Strothers'® and of Deutsch 
and Roberts" was not sufficient to resolve the 
two high energy gamma-rays. 

The relative abundance of the three gamma- 
rays can be obtained from the heights of the 
photoelectron lines if the effect of finite radiator 
thickness and the variation of the photoelectric 
cross section with energy are considered. The 
result is: 


Io.13 : I1.81 ° Io.845 = 20 : 30 : 100. 


COINCIDENCE MEASUREMENTS OF Mn‘** 


Coincidences between beta- and gamma-rays 
were observed with varying amounts of absorber 
between source and beta-ray counter as well as 
coincidences between gamma-rays. The beta-ray 
counter used was of the helium-filled bell type 
described in paper II. The gamma-ray counters, 
also helium filled, had platinum screen cathodes 
(paper III). The arrangement of the counters, 
designed to reduce the number of spurious 
coincidences due to scattered radiations to a 
minimum, has been described in paper V. The 
results are in good qualitative agreement with 
those of other observers* ® The number of beta- 


13 L. G. Elliott, M. Deutsch, and A. Roberts, Phys. Rev. 
63, 386L (1943). 

4 Since the photoelectric cross section in this energy 
range is not well known experimentally, the theoretical 
values of Hulme et al., Proc. Roy. Soc. 149A, 131 (1935), 
were used. 
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gamma coincidences per recorded beta-ray is 
independent of beta-ray energy above about 1 
Mev and has the value 1.00+0.05 X10- in the 
standard arrangement (paper V). When absorber 
is removed so that electrons of lower energy 
than 1 Mev are admitted into the beta-ray 
counter, the number of coincidences per recorded 
beta-ray increases rapidly and reaches a value 
of 1.92+0.0810-* when all absorber except 
the counter window (5 mg/cm?) is removed. 
The number of gamma-gamma coincidences per 
recorded gamma-ray is 0.92+0.07 X10-* in the 
standard counter arrangement. 


THE DISINTEGRATION SCHEME OF Mn‘ 


Since Mn** can decay only by negatron 
emission, the level scheme shown in Fig. 4a 
seems to be the only one consistent ‘with the 
experimental results. The fact that the high 
energy spectrum is accompanied by a single 
0.845-Mev gamma-ray is established by the 
excellent agreement of the observed coincidence 
rate of 1.0010-* and the efficiency curve for 
the counter used (paper V). If either of the two 
high energy gamma-rays found in the secondary 
electron spectrum were emitted in transitions 
leading directly to the ground state, the maxi- 
mum energy of the second beta-ray spectrum 
would have to be considerably higher than the 
value of about 1 Mev deduced both from the 
shape of the beta-ray spectrum and from the 
coincidence measurements. On the other hand 
there is not enough energy available for two 
high energy gamma-rays to be emitted succes- 
sively. The proposed disintegration scheme pre- 
dicts a value of 1.8810-* for the number of 
beta-gamma coincidences per beta-ray recorded 
without absorber, on the basis of the efficiency 
curve (Fig. 4, paper V) compared with the 
observed value of 1.92+0.0810-*. Similarly 
the predicted value of 0.99 X 10-* gamma-gamma 
coincidences per recorded gamma-ray is in 
excellent agreement with the observed value of 
0.92+0.07 X 10-*. Thus the disintegration scheme 
of Mn** may be considered well established. 


B. Co** 


In contrast to Mn** very little work has been 
done previously on the disintegration of Co**, 


Livingood and Seaborg'® have shown that the 
high energy positrons emitted with a half-life of 
70-80 days by cobalt obtained from deuteron 
bombardment of either iron or nickel must be 
due to Co**, since pure Co®* which was produced 
by alpha-particle bombardment of manganese 
emits much softer radiations. This was confirmed 
by producing Co** by fast neutron bombardment 
of nickel.'® 

Cook and McDaniel,’ using absorption and 
coincidence methods, reported that the positron 
spectrum is simple with a maximum energy of 
1.2 Mev, that the average gamma-ray energy is 
1.74 Mev, the energy of the highest energy 
gamma-rays 2.9 Mev, that the positrons are 
accompanied by gamma-rays, and that orbital 
electron capture seems to occur. 

In a Letter to the Editor® we have described 
a partial disintegration scheme for Co**. We 
shall now discuss the experiments leading to the 
adoption of that scheme and describe some 
other features of the disintegration of Co**. 


PREPARATION OF SOURCES OF Co** 


A study of the stable nuclear species in this 
region of the periodic table shows that it is 
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Fic. 5. Positron spectrum of radioactive cobalt separated 
from iron bombarded with deuterons. The curves represent 
spectra of three different sources. 


impossible to produce Co* free from other cobalt 
activities by any method used at present. It 
can be prepared by deuteron bombardment of 


16]. J. Livingood and G. T. Seaborg, Phys. Rev. 60, 
913L (1941). 

LL. G. Elliott and M. Deutsch, Phys. Rev. 63, 219A 
(1943). (See also paper VII, to appear soon.) 

17C. S. Cook and P. W. McDaniel, Phys. Rev. 62, 412 
(1942). 
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Fic. 6. Fermi plot of the positron spectrum of cobalt 
from (Fe+d). The plot represents the data shown in the 
top curve of Fig. 5. The low energy spectrum is due to Co**. 


iron in which case it is accompanied by Co, 
Co”, and Co*. Or it may be prepared by 
deuteron bombardment of nickel in which case 
Co®® and Co® are also produced. We have 
studied sources prepared by both of these 
methods. Most of the experiments were carried 
out with samples obtained from iron. The iron 
targets were treated as described in paper IV. 
The filtrate of the pyridine precipitation was 
freed repeatedly of elements of the hydrogen 
sulfide group (by using inactive copper carrier) 
and cobalt was precipitated as potassium cobalti- 
nitrite, and finally as oxide. Carrier material 
was added for all likely impurities, and all 
separations were repeated with additional carrier 
until all of the activity remained in the cobalt 
fraction. 

Sources for beta-ray studies were prepared 
either by deposition as sulfide as in the case of 
manganese or by electroplating on thin brass or 
copper foil. Gamma-ray sources were prepared 
by sealing the oxide into a suitable brass capsule. 

Samples prepared by deuteron bombardment 
of nickel were treated similarly. 


THE BETA-RAY SPECTRUM OF Co** 


Figure 5 shows the positron spectra emitted 
by sources of cobalt prepared by deuteron 
bombardment of iron. These curves do not 
represent the decay of a single source but 
spectra of sources prepared at different times 
before the experiments. It is seen that the spec- 
trum is complex and that the soft component 
decays somewhat more rapidly than the high 
energy radiation. Since the specific activity of 
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the oldest source was rather low, the corre- 
sponding positron spectrum is probably distorted 
by scattering in the source material. An estimate 
of this decay showed that very few, if any, of 
the soft positrons of Co*’ reported by Livingood 
and Seaborg!® were present. On the other hand, 
the conversion lines due to the gamma-rays of 
Co*’ with energies 0.119 and 0.131 Mev'!*!* were 
very pronounced, particularly in the older 
sources, and the intensity of these soft gamma- 
rays exceeded that of the gamma-rays of Co** 
and Co*®* many-fold. A discussion of Co*’ will 
appear in another paper. Figure 6 shows the 
positron spectrum of the strongest source as a 
Fermi plot. The plot is a straight line above 
0.48+0.01 Mev and indicates a maximum energy 
of the positrons of 1.50+0.05 Mev. After sub- 
tracting the high energy group, the low energy 
portion also appears to be represented by a 
straight Fermi plot. The maximum energy of 
0.48 Mev suggests strongly that it is due to the 
positrons'® of Co. This is supported by the 
more rapid decay of this group. Because of the 
presence of these positrons of Co’, the low 
energy portion of the Co** spectrum cannot be 
studied thoroughly, but it seems reasonably 
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Fic. 7. Secondary electron spectrum due to gamma-rays 
of cobalt from (Fe+d), obtained with moderate resolution. 
The rapid rise near Hp = 1000 is due to gamma-rays of Co*’. 


certain that any low energy group of Co", if it 
exists at all, must be much less abundant than 
the main group of maximum energy 1.50 Mev. 
No decay curves were taken of any of the sources, 
but the variation of the relative intensity of the 
two spectra (and of the gamma-rays as discussed 
in the next section) in sources prepared at 


18 E. H. Plesset, Phys. Rev. 62, 181 (1942). 
19 M.*Deutsch, A. Roberts, and L. G. Elliott, Phys. Rev. 
61, 389A (1942). 
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different times from the same bombardment is 
consistent with the values of 85 days and 65 days 
for the periods of Co** and Co**, respectively. 


THE GAMMA-RAY SPECTRUM OF Co** 


In a paper read before the Physical Society, 
Deutsch, Roberts, and Elliott'® have reported 
gamma-ray energies of 0.119, 0.131, 0.510 
(annihilation radiation), 0.835, 1.25, 3.40 Mev 
and a doubtful gamma-ray of 1.75 Mev from a 
mixture of long-lived cobalt activities produced 
by deuteron bombardment of iron. Figure 7 
shows an example of the data on which that 
report was based. It was pointed out at that 
time that the photoelectron line due to the 
0.835-Mev gamma-ray appeared wider than 
would be expected from instrumental effects and 
also that there was reason to believe that other 
gamma-rays exist of energies between 1.25 and 
3.4 Mev. No trace of the 0.20- and 0.22-Mev 
gamma-rays reported by Plesset'* was found, 
and it seems that the photoelectron lines ascribed 
by him to these gamma-rays were due to L and 
M electrons of lead ejected by the two known 
gamma-rays of lower energy definitely assigned 
to Co? which were the only ones found with 
energy less than that of the annihilation radiation. 
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Fic. 8. Study of secondary electrons in the neighborhood 
of 0.75 Mev with high resolution. Comparison of curves 
II, III, and IV shows the identity of the energy of the 
0.845-Mev gamma-rays of Mn** and Co**. Curves III, 
IV, and V Foe the varying relative abundance of Co** in 
sources of mixed cobalt isotopes. Note separation of K 
and L photoelectron peaks in curves I, II, and V. 


All secondary electrons of energy greater than 
0.85 Mev decay at the same rate, as shown by 
the fact that the shape of the distribution curve 
was found to be independent of the age of the 
source. The relative abundance of the annihila- 
tion radiation and of the peak due to the 0.84- 
Mev radiation decreases with increasing age of 
the source while the intensity of the low energy 
radiation due to Co’ decays much more slowly 
than the rest of the radiations. 




















‘ 


Fic. 9. Secondary electron spectra due to the high 
energy gamma-rays of Co**, observed with good resolution. 
Numbers indicate the location of the peaks and high 
energy end points, respectively. Solid circles show the 
location of the peaks due to high energy gamma-rays of 
Mn**, obtained with the same resolution. 


Since a 0.805-Mev gamma-ray has been 
definitely assigned'* to Co** and since a 0.845- 
Mev gamma-ray was found in Mn** which 
decays to the same product nucleus (Fe®*) as 
Co**, the photoelectron group at 0.75 Mev, 
previously ascribed to a gamma-ray of 0.835 
Mev, was investigated in more detail. The 
spectrometer was adjusted for better resolution 
and an arrangement was provided by which 
sources could be changed without disturbing 
the radiator. This eliminated any possibility of 
slight changes in the calibration of the spec- 
trometer due to shifts in the effective source of 
photoelectrons. Figure 8 shows the result of 
this study. A thin bismuth radiator was used 
rather than lead to take advantage of the larger 
photoelectric cross section and to increase the 
separation between the K and L photoelectron 
groups. Curves I, II, and V of Fig. 8 show this 
separation and illustrate the resolution used. It 
will be seen that both of the sources containing 
mixtures of Co®* and Co®* show two peaks. One 
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TABLE I, Gamma-rays accompanying the decay of Co*®. 
Relative intensities referred to that of the 0.845-Mev 
gamma-ray. If all disintegrations involve this gamma-ray, 
relative intensities also represent number of quanta per 
disintegration. Note that intensity of annihilation radiation 
is only 0.5. If all transitions were by positron emission, 
this intensity should be 2. 








0.510 
E (Mev) annihil. rad. 0.845 1.26 1.74 2.01 2.55 3.25 
Intensity 0.5 1.0 0.5 0.2 0.1 0.2 0.2 











of these corresponds accurately to the peak 
found with the source containing Co*® only. 
The other peak has an energy about 5 percent 
greater. The source prepared by deuteron bom- 
bardment of nickel shows a greater relative 
intensity of the Co®* gamma-ray, and a positron 
spectrum of the same sample showed a greater 
relative abundance of the Co*® spectrum. The 
higher energy peak obviously belongs to Co**, 
and it is seen by a comparison of curves II and 
IV in Fig. 8 that its energy agrees very closely 
(within 0.5 percent or better) with that of 
0.845-Mev gamma-ray of Mn**, 

Figure 9 shows a spectrum of secondary 
electrons due to the higher energy gamma-rays 
of Co**. Because of the large number of gamma- 
rays in this region and their low relative intensi- 
ties, it was necessary to take special precautions 
to reduce background effects (e.g., through 
coincidence counting) and to adjust the spec- 
trometer for very good resolution. In this manner 
it was possible to resolve five distinct gamma- 
rays where our previous measurements!® revealed 
only three. The upper curve shows the five 
groups of K photoelectrons, while the lower 
curve shows the Compton recoil groups due to 
the same five gamma-rays. The photoelectron 
group due to the highest energy gamma-ray 
appears somewhat wider than the others, sug- 
gesting the possibility of two unresolved peaks. 
Table I lists the energies of all of the gamma-rays 
of Co*®* together with their estimated relative 
intensities. 

Two of the gamma-ray energies, 1.74 Mev and 
2.01 Mev, are within a few percent of the values 
1.81 Mev and 2.13 Mev found in Mn**. However, 
they are certainly different transitions. The 
peaks due to the Mn** gamma-rays, obtained 
under the same geometrical conditions as those 
of Co**, are also shown in Fig. 9. There can be 





no doubt but that the difference between the 
energies is real. It should be clearly understood 
that, while we estimate the probable errors in 
the absolute energy values to be about two 
percent, this allows for possible errors in the 
calibration of the instrument and the effect of 
finite radiator thickness, while relative energy 
determinations within a narrow energy range 
should be reliable to about 0.5 percent. 


COINCIDENCE MEASUREMENTS OF Co** 


Coincidences between positrons and gamma- 
rays are studied preferably in the spectrometer 
since even with considerable precautions the 
annihilation quanta tend to introduce spurious 
effects when absorption methods are used. 
Coincidences between gamma-rays and positrons 
focused in the spectrometer were studied for 
several positron energies between 0.48 Mev and 
1.50 Mev. About three millimeters of lead had 
to be inserted between the source and the 
gamma-ray counter in order to absorb the very 
intense gamma-radiation of Co*’. In this arrange- 
ment the number of coincidences per recorded 
positron was found to be independent of positron 
energy and to have the value 4.0+0.210-. 
By measuring coincidences between beta- and 
gamma-rays of Mn**, Co’, and Co*®, whose 
disintegration schemes are known, it was found 
that the coincidence rate expected if each 
positron of Co** is accompanied by a 0.845- and 
a 1.26-Mev gamma-ray in cascade is 4.1+0.2 
x<10-*. This is in good agreement with the 
observed value. However, nearly the same co- 
incidence rate would be expected if each positron 
were accompanied by a single gamma-ray of 
about 2-Mev energy. Therefore the average 
energy of the gamma-rays accompanying the 
positrons was determined by observing the 
variation of coincidence rate as lead absorber is 
placed between source and gamma-ray counter. 
Figure 10 shows the result of this experiment, 
corrected for chance coincidences. On the same 
figure is shown the absorption (in the same 
geometry) of the gamma-rays from Co, whose 
disintegration is known to be accompanied by 
two gamma-rays of 1.1- and 1.3-Mev energy, 
respectively. The other dotted lines show the 
absorption expected with 2.1-Mev and with 
0.85-Mev gamma-radiation. The absorption co- 
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efficients found in this geometry are slightly 
larger than the accepted values because of the 
obliquity of the gamma-rays traversing the 
absorber. It appears that the observed absorption 
for the case of Co®® agrees very well with that 
expected for a 0.845-Mev and a 1.26-Mev 
gamma-ray in cascade, but would be quite 
inconsistent with the assumption of a 2-Mev 
gamma-ray. 

X-ray-gamma and gamma-gamma_ coinci- 
dences were also observed, but because of the 
unavoidable presence of other cobalt activities, 
particularly Co’, the results have little quanti- 
tative significance. 


THE DISINTEGRATION SCHEME OF Co** 


In the preceding section it was shown that the 
positron spectrum of Co** is accompanied by a 
1.26- and a 0.845-Mev gamma-ray as shown in 
Fig. 4b. The order of emission of the two gamma- 
rays is not directly established by experiment 
but the arrangement chosen in Fig. 4b seems 
very probable since a state of 0.845-Mev excita- 
tion energy in the Fe** nucleus was previously 
found in the disintegration of Mn** (Fig. 4a). 

The positron spectrum is certainly simple 
above 0.48 Mev and if any branching occurs at 
lower energies, the low energy spectrum is almost 
certainly not very abundant because of the 
competing orbital electron capture and the 
relatively low intensity of the high energy 
gamma-rays. It appears that the high energy 
gamma-rays accompany orbital electron capture 
almost exclusively. 

Any attempt to arrange these high energy 
gamma-rays into a level scheme is necessarily 
largely conjectural at the present time. A few 
features of the scheme can, however, be estab- 
lished with a certain degree of probability: It 
seems very unlikely that either the 1.74-Mev 
or the 2.01-Mev gamma-ray is emitted in a 
transition directly to the ground state from a 
level to which orbital electron capture takes 
place, since positron transitions to such a level 
would almost certainly be observed. The same 
argument holds—with less force—for the 2.55- 
Mev gamma-ray and for transitions of a 1.74- 
Mev gamma-ray to the 0.845-Mev level. The 
fact that the 0.845-Mev gamma-ray is about 
twice as abundant as the 1.26-Mev gamma-ray 





must have one of three possible explanations: 
1. Orbital electron capture directly to the 0.845- 
Mev state. 2. Other transitions of the same 
energy but between two excited states. 3. Other 
gamma-ray transitions, besides the 1.26-Mev 
gamma-ray, leading to the 0.845-Mev level. 

The first possibility seems extremely unlikely. 
In view of the great energy available for this 
transition (3.3 Mev), positron emission would 
certainly compete to a considerable extent with 
electron capture; and positron emission leading 
directly to this level is not observed. The second 
cannot be entirely excluded but seems somewhat 
improbable. The most reasonable assumption is 
the third. In Fig. 4c all of the states of Fe*®* are 
definitely established and the transitions be- 
tween them are indicated by solid lines, while 
the dotted lines indicate an assignment of the 
remaining gamma-rays to terms which are con- 
sistent with the above considerations. This 
particular arrangement is rather arbitrary and 
must necessarily remain so until further improve- 
ments in technique permit, for example, the 
study of coincidences between selected gamma- 
rays. 

The fact that the annihilation radiation is 
only half as abundant as the 0.845-Mev gamma- 
ray indicates that the total probability of orbital 
electron capture to all states is at least three 
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times as great as that of positron emission 
(except on assumption 2 above). The relative 
probabilities of the two modes of transition to 
the 2.11-Mev level alone cannot be determined, 
however, since part of the disintegration scheme 
is not yet fully understood. For the particular 
scheme indicated by dotted lines in Fig. 4c the 
probability of orbital electron capture leading 
to the 2.11-Mev level is about one-fifth of that 
of positron emission. 


C. DISCUSSION 


The existence of four excited states in Fe*® 
with excitation energies of 0.845+0.015 Mev, 
2.11+0.04 Mev, 2.66+0.05 Mev, 2.98+0.06 
Mev has been definitely demonstrated. Of these 
levels only the lowest one is excited both by the 
decay of Mn** and of Co**®. It can be seen, then, 
that neither disintegration involves all of the 
excited states which are energetically possible, 
and it is reasonable to assume that there are 
other states which are not excited in either 
disintegration. On the other hand the two sets 
of levels are not entirely unrelated since they 
have the lowest excited state in common. 

Inspection of the eight gamma-ray energies 
assigned to the Fe®® nucleus shows that there 
are several cases in which two of them may be 
added to obtain the energy of a third observed 
gamma-ray, or where one gamma-ray has just 
twice the energy of another. It has been shown 
in other disintegrations, e.g.,"%> Na™* and Br® 
(paper II) that such energy relationships do not 
necessarily mean that the same states are in- 
volved in the transitions. In Fe®®, too, we find, 
for example, that the excitation energy of the 
state to which positron emission takes place 
agrees, within the probable error, with the 
energy of the 2.13-Mev gamma-ray emitted 
following negatron emission by Mn**. 

Guggenheimer,” extending a suggestion of 
Rutherford and Ellis,” has pointed out that this 
behavior may not be fortuitous, since the energies 
of the gamma-rays emitted by some radioactive 
substances seem to occur in simple ratios. In 
Table II we list all of the gamma-rays emitted 


2°K. H. Guggenheimer, Proc. Roy. Soc. A181, 169 


(1942). 
Lord Rutherford and C. D. Ellis, Proc. Roy. Soc. 


A132, 667 (1931). 








by the Fe®® nucleus. Inspection of the column 
headed ‘‘E/0.425” shows that five of the eight 
gamma-ray energies are integer multiples of 
0.425 Mev, within the experimental uncertainty. 
It is particularly noteworthy that all integers 
from two to six are represented. The column 
headed ‘“Terms” shows the energies of the initial 
and final states of the transitions wherever 
these are known. 

The fact that the Fermi plots of the various 
beta-ray spectra are all straight lines within the 
experimental accuracy is rather surprising, since 
the positron spectrum of Co®* and probably also 
the high energy negatron spectrum of Mn** must 
be emitted in “forbidden transitions,’’ if this 
expression retains any significance with respect 
to a “Sargent diagram.” 

This means that all of the transition must 
involve differences of 0 or 1 in nuclear angular 
momentum between initial and final state (in 
the tensor theory), with or without parity 
change, the dominant matrix element being 
Sa or Sf1, respectively, in the notation of 
Konopinski and Uhlenbeck.'" The shape of the 
lowest (and possibly also the intermediate) 
spectrum of Mn** is somewhat uncertain, though. 
It is interesting to note that the ratio of the 
probability of orbital electron capture to that of 
positron emission predicted by theory” for an 
allowed transition in the case of Co*® is A/Aa 
=(.25 in close agreement with the value derived 
from the scheme shown in Fig. 4c. If the transi- 
tion is-first forbidden, this value could still be 
consistent with the theory for the tensor inter- 
action in certain transitions if the dominant 
matrix element is fa.” If the transition is of a 
different type, agreement with theory could be 
obtained only by revising the uncertain parts of 
the level scheme. 


TABLE II. Gamma-rays emitted by the nucleus Fe**, 








E Disintegration Terms E/0.425 
0.845 Mn*, Co 0.845-0 1.99 (2.0) 
1.26 Co 2.11 —0.845 2.97 (3.0) 
1.74 Co** ? 4.09 (4.0) 
1.81 Mn* 2.66 —0.845 4.26 - 
2.01 Co ? 4.73 — 
2.85 Mn** 2.98 —0.845 5.01 (5.0) 
2.55 Co** ? 6.00 (6.0) 
3.25 Co** ? 7.65 — 


| 
| 
| 
| 
| 





2 R. E. Marshak, Phys. Rev. 61, 431 (1942). 
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TABLE III. Possible angular momentum quantum numbers 
and parity of states of Mn, Fe, and Co. 











Excitation Angular 
Nucleus energy Mev momentum h/2x Parity 
Fe 0 0 e 
0.845 1 0 
2.11 6 e 
2.66 3 0 
2.98 3 0 
Mn* (4.21) 2 0 
Co** (4.12) 6 0 





It is possible to explain the absence of cross- 
over transitions in the disintegration scheme 
(Fig. 4c) on the basis of the known selection 
rules for beta- and gamma-ray transitions, as 
was done in other papers of this series. Table III 
shows a set of such angular momentum and 
parity assignments (the parity of a state is 
indicated by e if it is even, by o if it is odd). 
Of course, all the parities could be reversed 
without changing the character of the transition. 
The reasoning leading to these assignments is 
similar to that in paper V. The values proposed 
are not particularly convincing, and they are 
given only to show that no further selection 
rules are required to explain the present data. 

From the total disintegration energies of the 


schemes in Figs. 4a and 4b we can calculate the 
mass differences between the neutral atoms 
involved. We find Mn**— Fe®*= 3.98+-0.05 k 10-* 
aMU Co**— Fe®*=4.96+0.06 x 10-* aMU. 

Similarly one can calculate the threshold for 
the reaction Fe®*(p, 2)Co** from the disintegra- 
tion energies, taking into account the recoil 
energy of the product nucleus and the mass 
difference between neutron and proton ("—p 
= 1.26+0.02*10-* aMU). The threshold energy 
is thus found to be 5.47+0.12 Mev. This high 
value makes Fe** an interesting nucleus for the 
study of inelastic proton scattering, since it 
allows the use of high proton energies for the 
formation of excited states up to 5 Mev without 
competition from neutron emission. The much 
lower threshold'* for the reaction Fe®*(p, n)Co®*® 
should eliminate interference by protons scat- 
tered inelastically by Fe**. 
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The A-limiting process introduced by Wentzel and Dirac for the interaction of an electron 
with an electromagnetic field is applied to the interaction of a heavy particle (nucleon) with 
the meson field. One obtains in this way a more precise and relativistically invariant form of 
the theory which, in its classical interpretation, avoids all divergences of point sources. It is 
shown that the pseudoscalar and vector theories with a coupling constant f of the dimension 
of length have in this form the character of a weak coupling theory with no stable isobars 
existing for uf<1. The higher approximations for the interactions between nucleons are inves- 
tigated for the example of a Rosenfeld-Mgller mixture with the results that they are finite as 
long as the nucleons are treated to be at rest and small in comparison with the f? approximation 
for distances of the nucleons from each other sufficiently larger than f. 




























1. INTRODUCTION 





HE previous form of the theory of inter- 

action of a meson field with heavy particles 
(protons and neutrons, which are both called also 
“‘nucleons’’) suffered from the deficiency of a 
fundamental character in that the range of 
validity of the perturbation theory with respect 
to the coupling constant of this interaction could 
not be stated without introducing new assump- 
tions about the size of the nucleons. Since the 
higher approximations are divergent for point 
sources, one tried first to assume extended sources 
with a certain finite radius a. It turned out, 
however, that for the pseudoscalar and vector 
theories which alone give the right spin de- 
pendence of the nuclear forces, the perturbation 
theory never holds for values of a and of the 
coupling constant f with the dimension of length! 
which are required by experience. On the con- 
trary, for sources small in comparison with the 
Compton wave-length u-! of the meson and for 
af, the different approximation of strong 
coupling has to be applied. These strong coupling 
theories, however, had consequences in contra- 
diction to experiment.? Because of the existence 
of stable isobars with higher values of spin and 
charge, the highly charged nuclei should be 
unstable ; moreover, the magnetic moment of the 
proton and the free neutron should be just equal 
apart from their sign and should nearly cancel 



























. 1 In the following we use the natural unit where one puts 
=c=1. 

?W. Pauli and S. Kusaka, Phys. Rev. 63, 400 (1943), 
where other literature is referred to. 
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in the deuteron in contradiction to the experi- 
mental facts. 

One had the impression that the results of the 
perturbation theory are actually valid in a wider 
range than would follow from the model of the 
extended source. In the present paper* it is 
shown that this is understandable and can be 
formulated more precisely with the help of the 
so-called “‘\-limiting process’’ which was applied 
in the case of the electromagnetic field by 
Wentzel and Dirac in order to avoid the classical 
singularities of the self-energy of a point electron‘ 
in a relativistically invariant way. In Section 2 
it is shown that this limiting process can be con- 
sidered as the result of a natural generalization 
of the model of an extended source with respect 
to the reality conditions of the field variables in 
the momentum space, and that the main physical 
content of this formalism, which has the further 
advantage of preserving the relativistic invari- 
ance of the final results,5 can be characterized as 
making the constants of inertia of the degrees of 
freedom of the nucleon’s spin and the isotopic 
spin equal to zero. These latter constants were 


3 A short abstract of the contents of this paper was given 
in Phys. Rev. 63, 221A (1943), No. 25. 

‘Compare the literature quoted in the author’s report 
in Rev. Mod. Phys. 15, 175 (1943). 

5 In the present paper we do not give an explicit proof 
for this relativistic invariance, which can be done for the 
Hamiltonians of the type used in Section 3 with the help 
of Dirac’s formalism with several time coordinates in a 
way similar to Dirac’s proof for the case of electrodynamics. 
Compare also J. M. Jauch, Phys. Rev. 63, 334 (1943) where 
the application of the A-limiting process to the problem of 
the magnetic moment of the nucleons is treated. 





THEORY OF THE MESON FIELD 


defined by Bhabha® in his classical theory of 
neutral vector mesons, and he proposed also to 
assume their values equal to zero for the actual 
particles of nature.’ The discussion of the clas- 
sical model shows that no stable excited states of 
the nucleons exist in this theory if uf«1, a con- 
dition which certainly is fulfilled in nature and 
which characterizes the theory here discussed as 
a weak coupling theory in a more precise form. 

The Hamiltonians for the interaction between 
nucleons and meson field given in Section 4 do 
not contain anything new and are only prepara- 
tion for the investigation of the interaction 
between different nucleons in higher approxima- 
tions treated in the concluding Section 5. We use 
here the important method of successive canon- 
ical transformations which is due to Stueckelberg 
and Patry,® and we are particularly interested in 
the special mixture of a pseudoscalar and a vector 
meson field introduced by Rosenfeld and Mgller® 
which seems to be in better agreement with the 
empirical form of the nuclear forces than other 
assumptions about the meson field. While Rosen- 
feld and M@ller themselves merely guessed that 
the higher approximations for the nuclear forces 
would be relatively small if the distance r be- 
tween the nucleons is larger than the coupling 
constant f mentioned above, Stueckelberg’® 
showed that the higher approximations diverge 
for all distances r if the model of a point source 
is used because these higher approximations 
contain terms proportional to a~, where a is the 
radius of the nucleon. This result is in agreement 
with the later investigations from the standpoint 
of the strong coupling" theories in which the 
stable isobars of the nucleon play an essential 
role. 

In Section 4 we resume Stueckelberg’s inves- 
tigation of the Rosenfeld-Mgller mixture in order 
to apply the A-limiting process to this problem. 


*H. J. Bhabha, Proc. Roy. Soc. A178, 314 (1941). 

7It is an interesting question whether the d-limiting 
process can be generalized relativistically invariant in such 
a way that this constant is given an arbitrary value dif- 
ferent from zero and that divergences of theory are still 
avoided. My attempts to find such a generalization have 
not been successful. 

*E. C. G. Stueckelberg and J. F. C. Patry, Helv. Phys. 
Acta 13, 167 (1940). 

*C. Moller and L. Rosenfeld, Kgl. Danske Vid. Sels. 
Math.-Fys. Med. 17, No. 8 (1940). 

1” E. C. G. Stueckelberg, Helv. Phys. Acta 13, 347 (1940). 

" R, Serber and S. M. Dancoff, Phys. Rev. 62, 85 (1942). 
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For this purpose the use of the momentum space 
instead of the ordinary space in the description 
of the field variables is convenient, and it also 
simplifies the calculation itself. The result is the 
re-establishment of Rosenfeld and Meller’s orig- 
inal condition r>f for the validity of the per- 
tyrbation theory at least as far as the order of 
magnitude is concerned. It is true that the 
numerical factors in the higher approximations 
tend to diminish this range of validity of the f? 
approximation, but it has to be remembered that 
the other point of view, which is used in the 
above calculations and according to which the 
nucleons are considered to be at rest, restricts 
the validity of the results to r>M-', where M— 
is the Compton wave-length of the proton, which 
is practically of the same order of magnitude 
as f. It is therefore questionable whether the 
exact form of the higher approximations with 
respect to f* has any significance, and we merely 
stress their order of magnitude and the circum- 
stance that the A-limiting process is sufficient to 
make them all convergent so long as the nucleons 
can be assumed to stay at rest. 


2. THE 4-LIMITING PROCESS AS A GENERALI- 
ZATION OF THE EXTENDED SOURCE MODEL 


Choosing as the simplest example the inter- 
action of neutral pseudoscalar mesons with a 
heavy particle (nucleon) at rest in the origin of 
the coordinate system, we start with the well- 
known Hamiltonian 


H=3f (e+(vo)+utehdx 
+(4n)¥ f U(x)e-vedx, (1) 


where the real quantities r(x), g(x) are canoni- 
cally conjugate pseudoscalar fields; natural units 
with h=c=1 are used, yu is the rest mass of the 
meson in these units, f the coupling constant with 
the dimension of length, $e the spin of the nu- 
cleon, and U(x) the source function normalized 
according to 


fuw@av=1, (2) 


If brackets mean the Poisson symbols in the 
classical interpretation, and the commutators 
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multiplied by the imaginary unit 7 in the 
quantum mechanical interpretation, we have in 
both cases 


Lo(x), o(x’) ]=[x(x), x(x’) ]=0, 
(x(x), o(x’)]=6(x—x’); 
[ oi, a; | =— 2ox, 


(3) 
— ; ; (4) 
1, 7, k cyclic permutation of 1, 2, 3, 


and for the time derivative F of every observable 
F holds the equation of motion 


F=[H, F}. (5) 


In quantum theory oj, o2, o3 are the well-known 
spin matrices, while in the classical theory we 
can assume @ to be a unit vector. 

For the purpose of our generalization it is con- 
venient to pass from the ordinary space to the 
momentum space with the help of the Fourier 
transformation 


e(x) =(2m)-3 f g(ke**dke, 
(6a) 


n(x) = (2n)-* f pike dk, 


a(t) =(2n)-* f o(x)e*ax, 
(6b) 


p(k) = (20)-! f n(x)ei*¥dx. 


The rule (3) for the brackets is then equivalent 
to 


La(k), a(e’) ]=[b(e), p(k’) ]=0, 
p(k), g(k’) ]=5(k—k’). 


We further introduce the Fourier-transformed 
function v(k) of U(x) according to 


(7) 


U(x) = (2n)-* f v(e)e™, 
(8) 


o(b) = f Uae ae. 


The different choice of the normalization factors 
in (8) is convenient because the condition (2) is 


then equivalent to the simple statement 
v(0) = 1. (9) 


The Hamiltonian written in the momentum 
space is given by 


=} f { p(k) p(—k) + Retg(k)q( —) }dk 


— fo(—Be-ka(®)ak, (10) 
where 


ko=+(k2+y2)!. (11) 


In order to fulfill the conditions that U(x), 
a(x), g(x) have to be real, the quantities p(k), 
q(k), v(k) have to fulfill the reality conditions 


q(—k)=q*(k), p(—k)=p*(k), 
v(—k) =v*(k), 


(12) 


where a star denotes the conjugate complex. It 
is, however, remarkable that the reality of the 
Hamiltonian (and also of the total momentum) 
holds already, if the weaker conditions, 


v(k)q(—k) =v*(—k)q*(k), 
v(k) p(k) =v*(—k)p*(—k), 
p(k) p(—k) = p*(k)p*(—k), (13) 
q(k)q( —k) =q*(k)g*(—), 
v(k)v( —k) =v*(k)v*(—k), 


are fulfilled. The latter conditions are simplified 
by using the new variables 


q(k) =q(k)v(—k), plk)=p(k)v(k) (14) 
which satisfy reality conditions analogous to (13), 
q(—k)=q*(k), p(—k) =p*(k) 
and by using the real quantity 
G(k) =G*(k) =v(k)v( —R). (14a) 


The Hamiltonian in the new variables is 


H=} J [G(k) 7 {B(e)b(—k) + berg (k)G(—k) }dk 


if . 
omy -kg(k)dk. (15 
+ fe q(k) (15) 
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The value of the brackets is 


[h(k), g(k’) ]=G(k)5(k—k’). (16) 


The new form of the Hamiltonian makes it 
evident that models with different source func- 
tions v(k) belonging to the same function G(k) 
are equivalent, and that only the latter function 
has a physical meaning. According to (9) one has 
always 

G(0) =1, (17) 


while the particular case G(k) =const.=1 cor- 
responds to a point source. 

Whereas according to the stronger reality con- 
ditions (13) of the usual model of an extended 
source, G(k) is necessarily positive, the renounce- 
ment of the reality of the field variables in the 
ordinary space has led us to the weaker condition 
that G(k) has to be real only. 

It has to be shown that the particular choice 


G(k) =cos (Aoko —2-k) (18) 


has the remarkable property that an originally 
Lorentz invariant theory remains so if Ao, % is 
simultaneously transformed as a four vector. 
Finally, one has to carry through the so-called 
“\-limiting process,’’ namely, (Ao, 2)—0, which 
eliminates the particular choice of the \ vector 
from the final results and which makes all di- 
vergences of the classical model disappear pro- 
vided that the four vector Ao, 2 is always time- 
like, that is, if 


ho? >’. (19) 


As long as one has only to consider a single 
coordinate system, it is permissible to put 1=0; 
hence 


G(k) =cos Aoko =cos {Ao(k?+y?)!}. (18a) 


The A-limiting process avoids the (classical) 
divergences of the point-source model without 
introducing in the final results a finite extension 
of the source and without destroying the rela- 
tivistic invariance of a theory. 

In earlier papers we defined for the model of 
an extended source the reciprocal of its radius 
a by 


, 


dxdx 





@ *= 


f U(x) U(x’) 


/x—x’| 
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which in the momentum space can be written 


1 dk 
qa =—— few : (20) 
2’ k 


9 
T 


Maintaining this definition of a~' also in the 

case of a non-positive G(k) without any longer 

connecting a with the dimension of the source, 

we find for the particular choice (20) for G(k)” 
2 ~ 

o=— f cos {Ao(R?+ yu?) }dk=pJi(Aou). (21) 


wT “9o 


It is remarkable that the limiting process applied 
to the quantity a gives zero: 


Lim a'=0; (21a) 
oO 
while for the original model of a finite source, 
the transition to a point source, namely, G(k)—1 
means always a~'—«. Some physical conse- 
quences of the result (21a) are given in the next 
section. 


3. APPLICATION TO THE CLASSICAL THEORY 
OF PROTON ISOBARS (FREE GYRATION) 


As is well known there exists in the classical 
interpretation of the theory, under certain con- 
ditions for the value of the coupling constant and 
for the source function, the possibility of a con- 
tinuum of energy values for a heavy particle 
surrounded by its own meson field due to the 
presence of a free gyration of the field and the 
spin or isotopic spin of the heavy particle. This 
type of periodic motion was interpreted by Op- 
penheimer and Schwinger" as corresponding to 
discrete energy values in the quantum-theoretical 
treatment in which the total charge and the 
total angular momentum of the system have 
discrete values which are in suitable units integers 
and half-odd integers, respectively. They proved 
also that in the case of a strong coupling the 
treatment according to the correspondence 
principle agrees with the more rigorous quantum- 
mechanical treatment in the case of the charged 
scalar and the neutral pseudoscalar theory." It 


2 Compare for the evaluation of the integral, J. M. 
Jauch, Phys. Rev. 63, 334 (1943), Appendix. 

18 J. R. Oppenheimer and J. Schwinger, Phys. Rev. 60, 
150 (1942). 

‘4 The quantum-mechanical treatment in the case of the 
charged and symmetrical pseudoscalar theory was given 
by S. M. Dancoff and W. Pauli, Phys. Rev. 62, 85 (1942). 
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is the purpose of the following computation to 
investigate the conditions for the existence of 
such excited states (isobars) of the heavy particle 
if the \-limiting process is applied. It will be 
sufficient for it to give the classical treatment of 
the neutral pseudoscalar theory for the results 
in this case are typical also for more general 
cases. 

We start with the equations of motion which 
one obtains according to the application of the 
general rule (5) to the Hamiltonian (15) if one 


inserts for F the quantities g(k), g(k) and the 
components o1, ¢2, 3 of the spin of the heavy 
particle using the bracket expressions (4) and 
(16). In this way one obtains 





ag(k)_ 5(—b), 
al 
ap(k) af 
hence 
= uh 29(k) pay (22) 
and 
eu -— — f [oXk ]o(k)dk. (23) 


We suppose a free gyration of the spin around 
the x3 axis, according to 


o1=(1—C*)tcoswt, o2=(1—C?)!sin wt, 
¢3™= c. 


(24) 


C being a constant between —1 and +1. Insert- 

ing this in (22) we obtain 

§(0) = Ly | Os con wt +-by sin wf) 
=— ————~(k, cos in 

. v2 k?+ yp? —w* ’ " wai 


(25) 








ete" 
In the case 
w <p? (26) 


the denominator in (25) is always different from 
zero and the field g(x) in the x space decreases 
to zero for large distances from the heavy par- 
ticle. In the other case w?>v? the field in large 
distances has the form of a spherical wave and 
describes free mesons not bounded by the heavy 
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particle. In the following we shall therefore 
assume that the condition (26) holds. 

By inserting the expression (25) for g(k) in 
the right side of (23), the integrals can be sim- 
plified under the assumption that G(k) is 
spherically symmetric; that means that it 
depends only on the absolute value k of the 
vector k. In this case we obtain the following 
relation connecting the frequency w with the 
values of C 








4 ? 1 1 
w= fc f a) —~——}kdk. (27) 
3r 0 P+ pw—w ky? 


Before we discuss this relation we write down the 
analogous expressions for the energy and for the 
angular momentum. The latter is given by 


== + faxvelicte 


or in the momentum space by 


ni —1F ag(k) 1 
L= [ (60) | ak xk f+ bo. (28) 


It is in an integral of motion for a spherically 
symmetrical G(k) a specialization which has 
already been mentioned. Inserting in the ex- 
pression (15) for the Hamiltonian and in (28) 
for the angular momentum, the value (24) for e, 


“(—B) 


for p(k) one obtains for 





- 0 
(25) for g(R) and 
the energy 


f? 
E= B= f cwa-c| — 


Tv 





+p? —w* 


1 2w” 
— -— -— jes (29) 
R+p? (k?+p?—w?)? 
where 
f? kt 
{y= - f G(k)———-dk (29a) 
wo R?+p* 


is the self-energy of the nucleon corresponding to 
the ordinary static solution (C=1, w=0). The 
corresponding expression of the angular mo- 
mentum, which in our problem has the direction 
of the x; axis (L;=Z2=0, L3;=L) is 

k4 


2 2 
nvr wi-c) f eal Te ——dk+4C. (30) 
3m , 


+? —w?)? 











al 





iS 
K- 
3) 


or 


to 
he 
10- 
on 
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The integrals occurring in (27), (29), (30) may 
be evaluated in the following way. We split the 
integrand into one part without denominator, 
containing the integrals 


2 ¢” 2. 
N=- f G(k)kdk, at=- f G(k)dk, (31) 
T Ho | 


and another part which converges for a point 
source G(k)—1. In the latter part we shall insert 
G(k) =1 which is allowed both for the \-limiting 
process and for an extended source [G(k) 
positive ] with dimensions small in comparison 
with yo'!, that is, wa<1. We shall restrict our- 
selves to these two cases. It has to be emphasized, 
however, that in the former case for which our 
formulas will also hold, wa is in no way small. 
Using 


2 ¢* dk 
- f ————— = (y?—w")-}, 
rJo +o? 


2 f dk ute?) 
- =}(u—o*), 
wd, (e+ u2—w?)? 





and its specializations for w=0, one finds with 
the help of the indicated method, 


2 co) kA Ts 
— | G(k)——dk=N-——+', 
rv» k?+-y? a 


kt pu —w 


2 os) 
= [ G&)—_—ak=--~ + tw, 
ro hk? +p? —w? a 





2 k4 1 3 
— f G(k) -dk =-——(p?—w*)!. 
rvo (k?+ yu? —w*)? a 2 


Equations (27), (29), (30) get in this way their 
final form 





2 w 

w= -pc|——[at- (e—-0"), (32) 
3 a 

E=Evt( —C*) 


4 \—+u- (tw)! 34a), (33) 
a 
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Eo= -—(v-“+4), (33a) 
6 a 


1 t s 1 
L=-fra(t—C)|-—~t—w") +-C. (34) 
3 a 2 2 





For periodic motions the energy of which 
depends only on the amount of the resulting 
angular momentum, the frequency w is given by 
the general mechanical relation 


w=dE/dL. (35) 


Introducing the auxiliary function of the fre- 
quency 


W=wL—(E-—Ep), (36) 
this relation is equivalent to 
dW/dw=L. (37) 


The latter form can be checked easily in our 
case. We obtain first 


9 


1 f? w* 
W=-wC+—(1 = ©)|—— w+ (taf, (38) 
2 6 a 


and we note that according to (32) this can also 
be written 


1 i-C w 
W=-a(C+—— =—, (38a) 
2 2C 4C 


From (38) one derives 
1 2 w* 

aw=ac-| op ——s+ «|| 
2 3 a 


1 1 . @ 
+de C+-fra(t—C)|-——e—w")' |]. 
2 3 a 2 





The factor of dC vanishes according to (32), 
whereas the factor of dw is equal to L according 
to (34); hence the relation (37) is proved. 

The terms proportional to 1/a in the condition 
(32) for the frequency can easily be interpreted. 
By inserting (25) in (24) it follows that the terms 
proportional to 1/a in the equation of motion 
for @ are 


é=—ifta—[e, @]+---. (39) 


On account of the linear connection between 
g(k) and @ [Eq. (22) ] Eq. (39) holds generally, 
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and not only for periodic motions as would 
appear from its derivation. Bhabha" in his clas- 
sical treatment of the neutral vector mesons 
called the expression on the right-hand side of 
(39) the mechanical spin inertia term. The 
result a~'=0 of the A-limiting process means 
therefore that this process makes the spin 
inertia disappear; hence this process is in ac- 
cordance with Bhabha’s conjecture that for the 
real particles of nature the inertia terms will 
always be zero. 

We add the remark that for G(k) =cos Xoko, 
not only a~ does tend to zero together with Ao 
but also N defined in (31), and according to (33a) 
the self-energy of the nucleon then becomes 
Eo= —#f*x® which is small in comparison with 
the proton mass for the actual value of the 
coupling constant (fu)? ~ yo. 

We now turn to the discussion of whether Eq. 
(32) for w has a solution for —1=C=1 and w<u. 
We carry through the discussion for two cases. 
First we assume a small extended source cor- 
responding to a positive G(R) and to wa<1. In 
this case it is sufficient to retain the terms pro- 
portional to a~', and we obtain for the frequency 


3 a 
— 


2pC 





which can be fulfilled in the intervals in question 
for 


(fu)? > 3 (ua). (40) 


With the same approximation we obtain from 
(33), (34), (38): 
W=}twL=E—Ep, L=1/2C, 
3a 3a 


w=—L, E-E,)= 
f 2p 


The condition w <y implies.an upper bound for L. 
This is the old result of Oppenheimer and 
Schwinger which in the case of a strong coupling 
(fu)*>>wa agrees with the quantum mechanical 
result if we put Z?=/(/+1) with a _ half-odd 
integer / which is not too large. 

The second case we want to discuss is the 
result of the A-limiting process with a~'=0. 


(41) 





L?. 


1 See reference 6. In our units Bhabha's constant K is 
connected with our a“ by K/J = }f*a". 
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Putting x=w/yu Eq. (32) has here the form 
1—(1—x*)! 
1 = ¥(fu)*C———,, 
x 
which has solutions for O=x=1 and —1=C=1 
only if 
x: 3(fu)?>1, (42) 
where y is the maximum of the function 
[1—(1—x?)!]x—! in the interval O=x=1 which 
is only slightly larger than 1. This is not fulfilled 
in nature, where (fyu)?~ 75. Therefore in the 
theory based on the \-limiting process there does 
not exist a free gyration corresponding to a 
stable excited state of the nucleon with the actual 
value of the coupling constant, a result stated 
in the introduction. 

The theory of the scattering of a free meson 
can be treated in a similar way with this classical 
model. We do not give, however, this calculation 
in this paper because it contains nothing new in 
comparison with the analogous computation of 
Bhabha.'* The quantum theory of the meson 
scattering, especially the comparison of the 
results of the theory based on the A-limiting 
process with the theory of Heitler and Peng," 
needs further investigation. 


4. THE HAMILTONIAN FOR A MIXED PSEUDO- 
SCALAR AND VECTOR FIELD IN INTER- 
ACTION WITH SEVERAL NUCLEONS 


We must now generalize our Hamiltonian by 
including the isotopic spin, assuming the presence 
of several heavy particles and considering a 
mixture of a pseudoscalar meson and a vector 
meson. In the following we describe the nucleons, 
not with the help of a wave field in ordinary 
space and second quantization, but with the con- 
figuration space. The coordinates of the nucleons 
are denoted by z4 where capital Roman indices 
enumerate the different nucleons and run from 
1 to N if N nucleons are present. As in the last 
section we introduce in the interaction energy 


16H. J. Bhabha, reference 6, p. 333 ff. One has to identify 
his constant gs*/J with our f? and his constant 8 with our 
a~, The A-limiting process makes 8 = 0. The total scattering 
coefficient of the meson by the nucleon is identical with 
Bhabha’s expression (82), p. 337, if one substitutes for 
sin? @ its average } and includes a factor 2 which is due to 
the difference between the pseudoscalar theory here con- 
sidered and the vector theory treated by Bhabha. 

17W. Heitler and H. W. Peng, Proc. Camb. Phil. Soc. 
38, 296 (1942); W. Heitler, Proc. Camb. Phil. Soc. 37, 291 
(1941). 
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of the nucleons with the meson field first a 
general source function U(x) which degenerates 
into a 6 function for a point source and which 
will later be specialized according to the d-limit- 
ing process which re-establishes the relativistic 
invariance of the theory. 

The total Hamiltonian of the system consists 
of three parts Hs, Hy, Hu which describe the 
pseudoscalar mesons, the vector mesons, and the 
nucleons, respectively. For the sake of sim- 
plicity we neglect the interaction of the longi- 
tudinal vector mesons with the nucleons, and 
retain only the interaction of the transverse 
vector mesons and the pseudoscalar mesons. If 
we denote by a, 8 the Dirac matrices and put as 
usual 

7; —1Ba; 5,;= —1a2a3, oaes 
(43) 


Y¥5 = —181Q2Q3,; 


and denote with rt. (a=1, 2, 3) the components 
of the isotopic spin, one has according to Kem- 
mer!® and Bhabha!’ for the “‘symmetric’’ theory 


Hs=} J Zalwet+ (wen)? +e o')dx 


+(4r)'fs} 30 U(x—24)T04 
A,a@ 
X (84° V Ga— 54a} dx 
+2rfs?} >> U(x—2,) 
A, Ba 


X U(x—2p)ta*ta®ys*ys8dx+, (44) 
1 
tty =4 [ L)ne+—(0-m6) 
a Le 
+(¥ X $a)? +76." pdx 


+(49r)'ify | >> U(x—24)ra4 
A,a@ 


X {s4- (VX ba) — 74+ ta} dx 
+2rfy?} >> U(x—24)U(x—2zp) 
A, B.a 


x Ta! Tq®B4B?(s4 -s*)dx, (45) 


Hy=> 


A 





1 @ 
a4 -— —-+ Mp4 ' (46) 
4 oz* 

‘8 N. Kemmer, Proc. Roy. Soc. A166, 127 (1938); we are 
interested in his cases (c) and (d) and in the interactions 
proportional to his constants f, and ga. 

'*H. J. Bhabha, Proc. Roy. Soc. A166, 501 (1938). 
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In the limiting case of a point source, the theory 
is relativistically invariant, ¢ being a pseudo- 
scalar and ¢ together with go=yV-x a four 
vector. 7 is canonically conjugate to ¢ and = 
canonically conjugate to ¢ according to the 
bracket expressions (a, 8=1, 2, 3; i, R=1, 2, 3) 


[wa(x), p(x’) ]=5as5(x—x’), 


(47) 
[wai(x), gar(x’) ] = ba85%6(x —x’). 


Moreover, the components of the isotopic spin 
have the same commutation rules as the com- 
ponents of the ordinary spin, namely, 


[r1, T2]= —273, ---, (48) 


and the brackets of all quantities belonging to 
different nucleons are zero. With the help of 
these relations one obtains the equation of 
motion of all quantities by application of the 
general rule (5). 

As was shown by Kemmer, the terms of the 
Hamiltonian quadratic in the coupling constants 
are not unique, but there exist the two possi- 
bilities to choose for them 


2rf s? 7 U(x- Za) U(x— Zp) Ta* Ta? y54y58dx 


A, B,a 

or 

2rf s? Y U(x—24)U(x—2p) ta47a®(84-8”)dx 
A, B,a 


for the pseudoscalar meson, and 


2rfy? > Ul(x—24) U(x—2p) a4 ta? (y4 + yy?) dx 


A, B.a@ 

or 

2rfv?} dS U(x—24)U(x—2) 
A, B.a 


4 Ta? Tq®B4B7(s4 " s*)dx 


for the vector meson. The differences between 
these two alternatives are relativistic invariants 
in the limit of point sources.*° There is no a priori 
reason for one possibility or the other because 
both can be derived from a Lagrangian without 
explicit terms of the second order in the coupling 
constant with suitable independent variables. 


2? With the help of quantized wave functions ¥(z) and 
yt =y*B, pw, v=l, ---, 4, a4 = 1X0, these invariants can be 
written 2, [y* ivy F and 2 2 yt ivyuyw }, respectively. 
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The particular choice in the expressions (44) and 


(45) of these terms has only the a posteriori reason 
that it gives simpler results for the potential 
energy between two slowly moving nucleons, 
which in this case does not contain terms of the 
type 6(z4—2g). This situation, which has been 
noticed by different authors, does not seem 
entirely satisfactory. 

The form of the Hamiltonians Hs and Hy as 
integrals over the momentum space is given by 


Hs=} f Fel Palk)Pa(—k) +hetqa(k)ge(—k)} dk 


$ dX »(—k) exp (tk-Z,) 
—— —_ *Z Tat 
7 2? exp (tk-Z, 


X {t(s*-k)ga(k) — 54 pa(k) }dk 
fs? ; 
ae f =. [ik- (Za —Zp) |G(k) 


X tat taPys4ys®, (49) 


Hy=4 {Z| pelb)-Pa(—¥) 


1 
+—(K- pa(b)) (Ie: Pa(—)) 
7 


+ ko? qa(k) . qa( —k) 


— (K- qe(k))(K- qa(—)) pdk 


+ { 50-0) exp 2.) 
—- v(—k) exp (tk-Z4) ra4 
mv2 J A,a , 

x {is [hex qu(b) ]—14 pa() }d 
fv’ 


+— | > exp [tk- (z4—z2) ]G(k) 
4r’? A, B,a 


X ta47q2B4B7(s4-s*)dk. (50) 


The functions pa(k), ga(k) and pa(k), qa(k) are 
here defined as in Eq. (6), and they satisfy the 
canonical relation 


[pa(k), ga(k’) ]=5ap5(k —k’), 
[Pai(k), Pas(k’) -” 595; j5(k —k’), 


and reality conditions analogous to (13). The 
functions v(k) and G(k) are again given by (8) 


(51) 


and (15). Of course it is possible to define new 
quantities analogous to Eq. (14) of the preceding 
section which make it evident that only G(R) has 
a physical meaning. This will not, however, be 
necessary for the following. 

The part of Hy which describes the free mesons 
can be simplified by the canonical substitution 


1 m 
‘= p——k(k- p)+—k(k-q), 
p'=p—__k(k-p) + k(k-q 


‘ , (52) 
‘= q——k(k- q) -—k(k-p), 
q=q ie q uk p 
with the inverse formulae 
ete p') — ee @’) 
=p — e aus axes ° . 
p re p re q 
(52a) 
F i(k ) ic ’) 
q=q ie q uk Pp). 
One obtains in this way 
Hv=4{ Xtpe'(t)-pa!(—) 
+o? qa (Rk) -qa’(—k) }dk 
+f S-nepa-s 
—-- v(—k) exp (tk-z 
av2 A,a@ 4 
X Ta“ is4-[k- qa’(k) ] 
1 
=~" [pw OD) 
7 
— “c(e-@'(@) | 
k? 
fv’ 
--- >< exp [tk- (z4—Zz) G(R) 
4r? A, B,a@ 
X Ta4 T8488 (s4-s®)\dk+---. (50’) 


In this paper we shall deal particularly with 
the nuclear forces derived from the so-called 
Rosenfeld-Mgller mixture, which consists of a 
pseudoscalar and a vector meson with equal 
coupling constants. For the sake of simplicity 
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we shall assume also that the rest masses of the 
two kinds of mesons are equal. This may not be 
so in reality, but the purpose of the calculations 
carried through in the next section is not to 
compute the exact form of the nuclear forces, 
but merely to show the conditions for which the 
higher approximations of the perturbation theory 
regarding the nuclear forces are relatively small. 
For this purpose it may be sufficient to choose a 
simpler model than the actual case. Further sim- 
plifications arise from the approximation where 
the mass of the nucleon is considered infinite in 
comparison with the meson mass, or in other 
words where the recoil energy of the nucleon in 
all intermediate states can be neglected. In this 
case the four-component wave function of the 
nucleon can be reduced to a two-component wave 
function in a representation of the Dirac matrices 
where @ is diagonal, the two smaller components 
being negligible. All matrices which anticommute 
with 8, as for instance 7; and y, can be neglected. 
The four-component spin matrix s can be re- 
placed by the two-component spin matrix @#, and 
8 can be put equal to unity. The part Hy of the 
Hamiltonian is then unnecessary, and putting in 
our case 


fs=fv=g/n, (53) 


where g is dimensionless, one has 
H=Hst+Hy= 1fx | Pa(k) Pal —k) 


+ pa’ (Rk) Pa’ (—k) +ho*[ga(k)Ga(—k) 
+ qa'(k) -qa’(—k) ]}dk 


g 
+—— | > 0(—) exp (tk-z4)r.4is4 
rv2y A,a 


- {Kga(k) +[kX qa’(k) ]}dk 


2 





+ > exp [ik-(z4—Zzp) | 


4" A, B,a 
X G(R) ra47 a2 (s4-8®)dk+---. (54) 


As was shown by Stueckelberg,”! the simple 
properties of this particular mixture appear more 





_" E. C.C. Stueckelberg, see reference 10; the transforma- 
tion is made there in the ordinary x space. 
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clearly if one applies the canonical substitution 


q” = (t/u) {{kXq’]+kg}, 
p”’ = —i(u/k*){[kxXp’]+kp}, 


55) 
q’’ = —(t/n)(k-q’), 
p”’ =i(u/k*)(K-p’), 
with the inversion 
q’ = —i(u/k*) {(kXq" ]+kq’’}, 
p’ = (i/u){[kxp’]+kp”}, 
(55a) 


g=1(u/k?)(k-q’’), 
= —(i/u)(K-p”’). 


Moreover, the new functions fulfill the same 
reality conditions as the old ones. This substi- 
tution is convenient for the reason that the 
interaction energy depends only on q” but not 
on g”’. Because of the relation 


Pa(k) pal —k) + Pa’ (Rk) Pa’ (—*) 
+ ke qa(k) gal ™ k) + qa’ (k) qa’ ( —k)] 


k2 
=—[pe!(k)-Pa”"(—k)+Pe!(k)Pa’(—B) ] 
be 


2.,2 


Ue”) Ga”"(—B) +96") ge"(—#)] 





+ 


the field described by the new scalars pa”, ga”’ 
does not interact with the nucleon and can be 
split off. We do not write it down any longer, and 
we omit the double prime again in the following 
final result: 


1 k? 
H=- —Pa(k)-pa(—k 
5 J Ef Pa) -pel—®) 


2,,2 





a 


9 
«- 


qa(k) ) qa( —k) ik 


g 
4+ > v(—k) exp (ik- 
: 2 vo ) exp (tk-z,) 


X ra4e4- qa(k)dk 


on fz exp [ik-(z4—z5)] 


4rd A, B.a 





X G(R) ta4ra®(e4-@*). (56) 
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5. THE NUCLEAR FORCES BETWEEN NUCLEONS 
AT REST IN HIGHER APPROXIMATIONS 
FOR A ROSENFELD-M@LLER MIXTURE 

OF MESONS 


The method of Stueckelberg to find the inter- 
action energy between nucleons consists of the 
use of successive canonical transformations which 
eliminate step by step the field variables. The 
neglect of the motion (recoil energies in the inter- 
mediate states) of the nucleons is permitted as 
long as the distance between the nucleons is 
larger than their Compton wave-length M-!. 

Writing 


H=Hy+giit+¢*Hz (57) 
to express the dependence of the Hamiltonian on 


the coupling constant, we perform first a canoni- 
cal transformation 


H'=e"He* =H+[W, H] 
1 
sa hal [W,H]]+---, (58) 


where #7’ is the Hamiltonian in the new variables. 
While the middle term e'”He-*” holds only in 
quantum mechanics, the last form holds both in 
classical and quantum mechanics. We first 
assume W proportional to g according to W=gW, 
and arrange the result again in a power series of g 


H’=Hot+e{(W, HoJ+m} 


1 

+e{—[M, CM, HoJJ+lMm, H}+H| 
1 

+#\—[M, CW, (W:, Hol] 


1 
+ (Wi, (Wi, H+ (Ms, He] 1+ sss. (59) 


We choose W, in such a way that HA, is just 
canceled, namely, 


(Wi, Hy |j=-M, (60) 


and obtain 
H’=Ho+¢(4[ Wi, Hi] +H} 
+2°(30Wi, (Wi, Mi JJ+lWi, Hejj+---. (61) 


For the Hamiltonian (56) of the Rosenfeld- 


PAULI 


Meller mixture the condition (60) is fulfilled for 


1 
W,=-— v(k) exp (—7k-z,) 
awv2 ~ : 

ke? 


x Ta‘o4 - pa(k)dk. 





oe 
The term 3[ Wi, 1] which occurs in H,’ consists 
of two parts of which one is independent of the 
field variables while the other is bilinear in this 
variable. One finds 


4x? A, B, a, B, i, i 


Pi g 
—[W, WM Jj= -—-— a Ta Ta? at -@® 
2 


9 


x fex [ik-(z.4 ~2z5)]G(k)——dk 
ke? 


A 
(k’)v( —k) exp [i(k —k’) +24 
a ald fee )v(—k) exp [4( )+Za] 


12 


x ——pas(l’)4n( bald} (62) 
hou? 


| is an abbreviation for 


Here the coefficient : 
at, B7 


A 
[A Patriot 
at, B7 


A A 
= —20%jbap—26ijtes (63) 


(with the notation r12= —721=73, °° * } O12= — 02 
=0o3, -::). Because of the particular choice of 
H. a certain cancellation occurs between the last 
term in (62) and Hz since 

1 k? 1 


wu? Rot? eo? 
with the consequence that in the limit of point 
sources, G(k)—1, no interaction of the type 
§(24—Zg) appears in the result. We get finally 


for the part H3,o which is independent of the 


field variables—generally in the symbol Has 
the first index is equal to the power of g with 
which it is multiplied, and the second index 
denotes the degree in the field variables—the 
result 
at A tase f 
2H 9 3 ae Tao -@ J pie 
Xexp [ik-(z4 —Za) |dk. 
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orders of the factors 


THEORY OF 


For A ¥B it is always allowable to put G(k) = 1. 
If we use 


1 exp [ik- (za. ~2)1 


Beccs. (—urap) 
— ———, (64) 
2r* k? +p? TAB 





rap=|Za—Ze| being the distance between the 
two nucleons and take into account that every 
pair A, B occurs twice in the sum }° 4,2 we find 
for the potential between two different nucleons 
A, B as far as it is proportional to g?, the well- 
known result 


Vap=2 Doe Ta*Ta8o4 -o*® [exp (—yuraz)/rap }. (65) 


For A =B we obtain the self-energy the value of 
which in the limit G(k)—1 is given by 


9 1 
Vaa=—¢(--«), (66) 
2 \a 


if, according to quantum theory, )-« Ta? = (4)? = 
is inserted. For the case of the A-limiting process 
we have to put a~'=0; hence 


9 
Vaa=— oe (66a) 


In the discussion of terms of higher order in g 
we are mostly interested in the terms linear in 
the field variables, the lowest order of which is the 
third. These terms are according to (61), (62) 
given by 


3 AB 
i ee, 2 ‘ ‘| 


(2m)! 4 pail at 


xf “ exp [ik- (Z1—Zp) \dk 


2 


xf exp (—ik’-z,4) ——0(k") pai(k’) dk’, 
ky “Kh 


AB A 
“|=z| rots (67) 
at 8,j Lat, B7 


In the case A = B one has to hermitize the right 


where 


side by taking the arithmetical mean of the two 


A A, .A > 
od . and rg4¢;*. One finds 


THE MESON FIELD 


in this way for A=B 


AA 
° |=o, (67a) 


at 


for A#B 
AB 
. ]=2t extort extort) rat 
at 
+o;8 (734 rT," 


where (i, 7, k) and (a, 8, y) are cyclic permuta- 
tions of (1, 2, 3). 

These new linear terms in the field one has 
now to transform away by a new canonical trans- 
formation of the type (58) with the help of a 
W=g°W; which satisfies the condition 


(W;, Ho)=—H31 


—t,47g")} (67b) 


and is given by 
1 


W;= = 
(2x)! Ss 


fe 


x [exp (ik’ -Z4) 
e ko’ 





“9 


al 





exp [ik- (Za —Zp) \dk 


2 


0(—R')gai(k’)dk’. (68) 
22 





This transformation gives rise to a part H¢ o of 
the new Hamiltonian which is contained in 


HH, * =4[ W;, H;, 1}. One finds” 


ABYCD 
oH, o= 3¢° } ; | wi |Feafanfeo (69) 


A, B,C,D,a,i @&t 


with the definitions 








1 G(R) 
fas= -f be exp [tk-(zZ4—Zz) |dk, (70) 


3 


ee = fow— 
ko! we 


Xexp [tk-(z4—Zaz) Jdk. (71) 


For the limiting case of a point source the 
integral f4g was evaluated above as 


fan=exp (—uraz)/raz for A¥B, (72) 
faa=a'—up. 


* Compare E. C. C. Stueckelberg, reference 10, formulas 
(2.10) and (2.11). 
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In a similar way one finds 


1 2 
Fan=—(1- ) exp (—uran for A+B, 
2u MTAB 


1 2 
Faa -—(s-—). 
2u ua 


The result (69) means in general a three- and 
four-body interaction besides the two-body inter- 
action. In the case of only two particles present 
the expression (69) gives an additional potential 
energy given by 





(73) 








7) g® 3 exp (—2ura B) 
Vap=—| Cas ; 
# TAB 
2 exp (—3yuraz) 
+dan(1 - ) - | (74) 
MrAB vB 


where A #B and 


ABYAB ABYBA 
can=Z| “| { iun= | ‘| | (75) 
a,i at at a,t at at 


and where we put again a~'=0 according to the 
\-limiting process. 

The evaluation of c4g and dag with the help 
of the expression (67) gives, with convenient 
abbreviations, S4s=(@4-¢8), TAB= doa Ta*Ta?: 


Cap=8{3(3—Sap)+3(3—Taz)t+4sanras}, (76) 
dap= —8{(3—Saz)Tap 


+5s4e(3—7az)—4S8aptap}. (77) 


22 E.C. C. Stueckelberg, see reference 10, Eq. (3.12) gives 
in both cagand daganadditional term —64(3—s48)(3—raz) 
which I have been unable to confirm and which seems to 
be erroneous. This additional term causes also an over- 


: ° ° it 
estimation of the numerical value of V4g by Stueckelberg. 
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For the ground state of the deuteron one has 
Tap= —3, Sap=1, Cap=3.32, dap= —3.32, and 
V"/V small for 3(2g)*/(ur)?<«1 as was indicated 
in the introduction. 

While in the classical theory the only terms of 
the final Hamiltonian which give rise to an inter- 
action between nucleons are those which are 
independent of the meson field, it is different in 
quantum theory. In this theory it is not possible 
to put simultaneously p(k) and ga(k) equal to 
zero, and terms in the final Hamiltonian which 
are either bilinear in the ~.(&) or bilinear in the 
ga(k) give in general a finite contribution to the 
nucleon interaction due to the zero-point energy 
of the field oscillators. There exist such terms of 
the order g‘ which are generated by a g?W2 which 


eliminates g?H:2,2 given by the first term of (62), 
and which are computed by different authors. It 
is interesting that this interaction energy of the 
order g* disappears again according to the new 
theory of Dirac, where mesons of negative energy 
are introduced in the intermediate states. This 
theory brings forth a greater similarity with the 
classical theory than the older form of the 
quantum theory, and particularly all effects due 
to the zero-point fluctuations of the field oscil- 
lators are canceled in the new theory. It would 
be possible to check this new theory if experi- 
mental tests for or against the fourth-order terms 
of the interaction energy were available. 

On the other hand this new hypothesis of 
negative-energy mesons was not necessary in 
order to make the theory convergent in the ap- 
proximation where the heavy particles are treated 
non-relativistically. In this approximation the 
convergence could be achieved with the \-limiting 
process alone. 
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It is pointed out that the cross sections for the production of mesons in collisions between 
heavy particles given by Hamilton, Heitler, and Peng are so large that a heavy particle crossing 
an atomic nucleus must be expected to collide several times inside the same nucleus. This 
hypothesis accounts for cloud-chamber observations showing the production of several mesons 
in one point. Anticoincidence experiments of Janossy and Rochester can also be interpreted 
easily by the theory of HHP when the occurrence of multiple collisions is assumed. 





I. INTRODUCTION 


N a recent publication Hamilton, Heitler, and 
Peng! (quoted as HHP) investigate theo- 
retically fast collisions between neutrons and 
protons. It is found that in fast collisions between 
heavy particles mesons are emitted. The cross 
section for the emission of mesons is found to be 
of the right order to account for most cosmic-ray 
phenomena. In the following we attempt to add 
to one point of the theory of HHP. 

Cloud-chamber photographs of penetrating 
showers obtained by various observers? fre- 
quently show small groups of mesons diverging 
from a common center. These photographs are 
most easily interpreted on the assumption that 
mesons are produced simultaneously in small 
groups. 

A typical photograph of this kind obtained 
recently by Janossy, McCusker, and Rochester‘ 
is reproduced in Fig. 1. The photograph shows 
three tracks diverging from a point above the 
chamber. Each of the three tracks is passing 
through a lead plate, placed across the chamber, 
but none of the particles starts a cascade while 
passing through the plate. Thus the tracks cannot 
be assumed to be electrons. The simplest inter- 
pretation of the photograph is to assume that at 
least two of the three particles are mesons pro- 
duced simultaneously in one point. The third 
particle is either a third meson or the primary 
particle which created the group. 

The theory of HHP assumes that a fast heavy 


1 Hamilton, Heitler, and Peng, Phys. Rev. 64, 78 (1943). 
*Braddick and Hensby, Nature 144, 1012 (1939); 


G. Herzog and W. H. Bostick, Phys. Rev. 58, 218 (1940). 
3 W. M. Powell, Phys. Rev. 58, 474 (1940) ; 60, 413 (1941). 
‘W. H. Bostick, Phys. Rev. 61, 557 (1942); Janossy, 

McCusker, and Rochester, Nature 148, 660 (1941). 
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particle produces mesons one by one in succession. 
Multiple processes are mentioned as a possible 
source for the groups but are not dealt with by 
HHP. We wish to point out in the following that 
the groups of mesons observed can be accounted 
for in a simple way without considering multiple 
processes in the strict sense. 


Il. MULTIPLE COLLISIONS 


The cross section for a fast neutron or proton 
to produce a meson with the energy greater than 
E in a collision with another heavy particle is 
according to HHP 


@(E£) =1.5X10-*/E cm’. (1) 


All energies are measured in units of yuc*?, where u 
is the mass of the meson. 





Fic. 1. Meson shower, showing three penetrating particles 
passing through a lead plate. 
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_ For moderately large values of the energy E 
the cross section given in (1) is so large that a 
fast heavy particle passing through a nucleus is 
likely to collide several times during the passage 
through the nucleus; several mesons may there- 
fore be emitted in a collision with one nucleus. 
The collision of a fast heavy particle with a 
nucleus can be considered schematically as 
follows. 

The radius of a nucleus of the weight A can be 
assumed to be r4=0.53(e?/mc*)A* (compare 
Heisenberg®). The centers of the heavy particles 
forming the surface of the nucleus thus have an 
average distance 74 from the center of the nucleus 
given by (2), 


F4 =0.53(¢%/mc*)(A*—1). (2) 


Geometrically a collision between a constituent 
of the nucleus and the fast particle can take place 
when the fast particle passes at a distance less 
than x =74+ rg from the center of the nucleus, rz 
being the radius of a disk with the area (£) 
whence mrg?=(E). The heavy particle passing 
through the nucleus has a chance to pass between 
the nuclear constituents without a collision. The 
probability for a passage without collision can be 
estimated roughly as 


p=exp [—®(E)/(x74°/A) J, (3) 


and therefore the effective cross section for the 
emission of at least one meson of the energy £ in 
a collision with a nucleus of the weight A is 
estimated roughly as 


(A, E)=(1—p)a(re+Fa)?. (4) 
Inserting numerical values one finds 
@(A, E)=1.5X10-*/E 
X (1+aE!)?(1-—exp[—b/E]). (4a) 


The constants a and 6 are functions of A only: 
Some values are given in Table I. Note that 
(A =1, E)=4(£) and that for large values of 
E®(A, E)—A®(E) as it must be expected. 











TABLE I. 
Pb Al Oo H 
a 1.05 0.43 0.33 0 
b 178 140 145 2 








5 Heisenberg, Leipz. Ber. 89, 369 (1937). 
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Fic. 2. Mean free path of fast heavy particles in lead. 
Curve I, 1/Ry, packing neglected. Curve II, 1/Ry, com- 
pact nucleus. 


The reciprocal mean free path 1/R= N®(A, E) 
as function of E for heavy particles in lead is 
plotted in Fig. 2 with the full line. N is the 
number of atoms per cc. The broken line marked 
I in Fig. 2 represents the reciprocal mean free 
path, 1/R1= NA®(E), one would expect if the 
packing of the nuclear particles could be neg- 
lected. The broken line marked II represents 
1/Ru = Nx(re+fa)’?, the reciprocal mean free 
path one would expect if we neglect the possi- 
bility of a fast particle traversing a nucleus 
without encounter, i.e., regarding the nucleus as 
compact. 

It is seen from Fig. 2 that for the emission of 
mesons with energies up to 5X10° ev (E=50), 
the lead nucleus can be regarded as compact, 
while for the emission of mesons of energies above 
10" ev (E=1000), the nuclear particles can be 
regarded as independent. 

The average number m4 of mesons with energies 
exceeding E emitted by a sufficiently energetic 
fast particle while traversing a nucleus of weight 
A is expected to be 


n4=A(E)/*(A, E) 
=A(1+aE})-*(i+exp[—b/E])-. (5) 


The multiplicities thus obtained for lead are 
given in Table II. For lighter nuclei the multi- 
plicities are somewhat smaller. The values ob- 
tained for m4 include pseudoscalar and transverse 
mesons as well as neutrettos. 

The argument leading to Eqs. (4) and (5) 
assumes the nuclear particles to be independent. 
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This is of course not correct and in the accurate 
treatment interference would have to be taken 
into account. In the case of the emission of slow 
mesons the effect of interference may be serious. 
In the case of fast mesons, however, the de 
Broglie wave-lengths of the primary particle and 
of the emitted mesons are small compared with 
the average distances of the nuclear particles and 
the frequencies of both the primary particle and 
the emitted mesons are large compared with the 
nuclear frequencies. It is therefore likely that the 
results obtained on the assumption of independent 
nuclear particles are at least qualitatively correct 
for the emission of fast mesons. 

We conclude therefore that the theory of HHP 
leads us to expect that mesons are produced in 
small groups. 


III. COMPARISON WITH EXPERIMENTS 
A. Cloud Chamber 


The values of the multiplicities m4 given in 
Table II are, of course, only valid for primaries 
having sufficient energy to produce the required 
number of mesons. Primaries having less energy 
will be stopped completely inside the first nucleus 
with which they collide closely. The energy re- 
quired to pass through a nucleus is of the order 
of 200. The energy required to pass through the 
atmosphere is according to HHP of the order of 
10,000; the average energies of primaries which 
succeed in traversing the atmosphere must be 
expected to be of the same order. Thus it is ex- 
pected that most of the primary particles 
reaching sea level have sufficient energy to 
traverse a lead nucleus. 

The multiplicities m4 seem, however, to be too 
large to be compatible with the experimental 
evidence from cloud-chamber photographs. The 
predominance of small groups on photographs 
containing meson showers may be explained by 
assuming that most of the fast heavy particles 
near sea level are recoil particles produced by the 
primaries. The average energy of the recoil 
particles must be assumed to be small; and most 
of the recoil particles will have energies only 
sufficient to produce small groups. 

In Fig. 3 a photograph of a big shower obtained 
by McCusker is reproduced. As the chamber was 
controlled with a heavy bias. towards penetrating 
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Fic. 3. Large shower obtained with an arrangement giving 
strong bias towards penetrating showers. 


showers, it is not unlikely that this photograph 
represents a penetrating shower produced by a 
very energetic primary proton or neutron. The 
tracks seen on the photograph are, however, so 
numerous that it is impossible to decide whether 
or not some of the particles are mesons. 

Shower photographs of Blackett and Occialini® 
and of Anderson’ give evidence for the occurrence 
of slow heavy particles in showers. Janossy, 
McCusker, and Rochester* have suggested that 
the occurrence of slow heavy particles is con- 
nected with penetrating showers. We note that 
the occurrence of slow heavy particles in pene- 
trating showers is well accounted for by the 
theory of HHP in terms of slow recoil particles. 


B. Counter Experiments 
1. Transition Effect of Penetrating Showers*® 


Using a counter arrangement as shown in 
Fig. 4 we have recorded coincidences involving 
TaBLeE II. Number of mesons with energy greater than E 


emitted by a fast heavy particle crossing a lead nucleus. 


E 5 10 100 1000 
nA 18 11 1.7 1 














6 Blackett and Occialini, Proc. Roy. Soc. A139, 699 
(1933). 

7C. D. Anderson and S. H. Neddermeyer, Phys. Rev. 50, 
263 (1937). 

*L. Janossy and Ingleby, Nature 145, 511 (1940); L. 
Janossy, Proc. Roy. Soc. A179, 361 (1941). 
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Fic. 4. Arrangement for observing penetrating showers. 


the discharges of more than one counter out of 
each of the counter trays B, C, and D. The 
counter trays were separated by lead absorbers 
15 cm thick, and the whole arrangement was 
surrounded with an absorber 50 cm thick. It was 
shown that most of the coincidences were due to 
penetrating showers. The rate of penetrating 
showers was found to increase when a lead 
absorber T was placed close above the top tray B. 
The increase was considerable for T a few cm 
thick, while no increase was found for thicknesses 
larger than five cm of lead. More detailed 
measurements of the transition effect will be 
reported shortly by Rochester and Janossy. 


2. Production of Penetrating Showers by Non- 
Ionizing Radiation (Jénossy and Rochester®) 


The arrangement is shown in Fig. 5. The lower 
part of the arrangement consisted of a coinci- 
dence arrangement sensitive to penetrating 
showers. The absorber T was surrounded with a 
set of anticoincidence counters A. Anticoinci- 
dences A—BCD, that is, coincidences BCD not 
accompanied by the discharge of any of the 
counters A, were recorded. These anticoinci- 
dences were largely due to penetrating showers 
started by non-ionizing particles falling on 7. 

To obtain information about the mean free 
path of the non-ionizing radiation giving rise to 
penetrating showers an absorber = was placed 
over the anticoincidence counters. It was found 
that five cm of lead placed above the counters A 
did not noticeably affect the anticoincidence rate. 
However, a lead absorber 35 cm thick placed over 
the anticoincidence counters cut down the rate of 

* L. Janossy and G. D. Rochester, Nature 150, 633 (1942); 


Proc. Roy. Soc. A (in press) ; J. G. Wilson, Proc. Roy. Soc. 
A174, 73 (1940). 


anticoincidences considerably. It was thus con- 
cluded that the mean free path of the non- 
ionizing radiation well exceeded 5 cm of lead, 
Thus the non-ionizing radiation could not be 
assumed to consist of photons, and it was sug- 
gested that it might consist of neutrons. 

To avoid confusion we note that the term mean 
free path is used in the sense as the average 
distance a particle has to travel before encounter- 
ing a collision of some specified kind. This mean 
free path is in general different from the range 
of a particle. If a particle is capable of suffering 
several collisions before being brought to rest, the 
range exceeds the mean free path considerably. 


C. Interpretation of the Counter Experiments 


The counter experiments described above were 
interpreted by HHP in terms of incident protons 
and neutrons. While we agree with this interpre- 
tation, we wish to point out that considerable 
difficulties are removed by introducing our 
hypothesis that a heavy particle collides several 
times in one nucleus. 

1. If the collisions of a fast primary were dis- 
tributed at random along its path and not 
crowded together into nuclei, then the mean free 
path for the production of mesons of energies 
greater than E could be obtained from Eq. (1) as 


Ri =1/NA®(E) =E/10 cm Pb. (6) 


Thus a fast particle traversing a few mm of lead 
would be expected to give rise to secondaries 
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Fic. 5. Anticoincidence arrangement (JAnossy 
and Rochester). 
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which would emerge together with the primary. 
This, however, is clearly incompatible with the 
interpretation of the anticoincidence experiment 
in terms of neutrons described in Section III, B2. 
The anticoincidence experiment shows that the 
non-ionizing particles have a good chance of 
penetrating five cm of lead without an encounter 
leading to ionizing secondaries emerging out of 
the lead. 

Considering, however, the crowding of colli- 
sions, the mean free path of a primary has to be 
obtained by using Eq. (4). One finds 


R=(1/N)®(A, E). (7) 


Inserting numerical values one sees that R is 
much greater than R;. R, however, decreases 
rapidly with decreasing energy E, and the value 
Ry of R corresponding to the mean free path as 
observed in the anticoincidence experiment is 
obtained as follows. 

Collisions giving rise to mesons of energies so 
low that they are stopped inside the absorber do 
not prevent a neutron from giving rise to an 
anticoincidence.” Such collisions can therefore be 
neglected when considering anticoincidence ex- 
periments. Roughly one can neglect all collisions 
corresponding to values of E smaller than Ep 
when the value of E, is chosen in such a way that 
the mean free path R= R, obtained from (7) for 
E=E, is equal to the range of mesons of the 
energy Eo. For lead one finds 


Eo~2, Ro~7 cm Pb (8) 


in good agreement with the experiment. 

2. The transition effect of penetrating showers 
described in Section III, B1 can be interpreted in 
terms of protons and neutrons giving rise to 
penetrating showers in the absorber 7. Saturation 
is expected for thicknesses somewhat exceeding 
Ro. This is compatible with the experimental 
results. 

#0 A neutron, however, can be transformed into a proton 
while suffering a low energy collision, and thus it can be 
prevented from giving rise to an anticoincidence even if the 
mesons created in the collision remain inside the absorber. 
The loss of neutrons due to this process is, however, 


compensated by the gains of protons which are converted 
into neutrons in low energy collisions. 
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Comparing the results of the two experiments 
described in Section III, B, one finds that five cm 
of lead put into the path of the non-ionizing 
radiation has little effect on the anticoincidence 
rate, and therefore the mean free path of the 
neutral particles appears to be well above five cm 
of lead. Regarding the ordinary transition effect 
(Section III, B1), we note that most of the 
transition effect is over after five cm of lead. 
Thus it appears as if the ionizing radiation pro- 
ducing penetrating showers had a shorter mean 
free path than the non-ionizing radiation. If the 
two radiations consist of neutrons and protons as 
we think they do, this may indicate that the 
mean free path of protons is somewhat less than 
that of neutrons. 

On the other hand the difference between the 
two experiments can probably be accounted for 
without assuming different mean free paths for 
protons and neutrons in the following way. 

Penetrating showers are known to be ac- 
companied by extensive air showers. A neutron 
falling on the anticoincidence arrangement 
(Fig. 5) may fail to produce an anticoincidence 
when accompanied by a sufficiently dense air 
shower as the air shower is likely to discharge 
some of the anticoincidence counters. 

An absorber = five cm thick placed above the 
anticoincidence counters has a twofold effect. 
First, some of the neutrons which would be 
recorded without the absorber = give rise to 
mesons while passing through it and thus are not 
recorded. Secondly, some neutrons originally 
accompanied by air showers may be stripped of 
the shower while passing through >. Thus 
neutrons which are stripped of the accompanying 
cascade and happen to pass through 2 without 
encounter are recorded additionally because of 
the presence of 2. The two effects may compen- 
sate each other to some extent, and it is quite 
likely that this accounts for the small absorption 
observed in the first five cm of lead. 

I am indebted to Professor W. Heitler for 
valuable comment and for communicating his 
results, and also to Dr. G. D. Rochester for valu- 
able discussion. 
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Adiabatic Compressibilities of Some Aqueous Ionic Solutions and Their Variation 
with Indicated Liquid Structure of the Water* 


Victor B. Corey 
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The adiabatic compressibilities of aqueous solutions of twenty-seven strong electrolytes, 
each at a single concentration, were determined at 25°C from sound velocity measurements 
at a supersonic frequency (1015.24 kc/sec.). The change in compressibility per mole fraction 
of solute shows a satisfactory correlation with x-ray and other indications of water structure 
change. Adiabatic compressibility values themselves show a remarkable correlation with 
corresponding values of the partial molal volume of the solvent water. Evidence strongly 
supports the conclusion that water has a liquid structure which becomes more highly coordi- 
nated and compacted with the introduction of ions. 

















































INTRODUCTION 


HIS report is upon the adiabatic compressi- 

bility of pure water and of aqueous solu- 
tions of twenty-seven strong electrolytes, and an 
interpretation of the results in terms of the 
liquid structure of the solvent water. Isothermal 
compressibility measurements by Amagat! and 
Bridgman? showed a minimum in the compressi- 
bility of water, and a dependence of compressi- 
bility values and the temperature of the mini- 
mum upon pressure. Numerous piezometric 
measurements have been made on water and its 
ionic solutions.* By modern supersonic methods, 
the adiabatic compressibilities of many organic 
liquids have been obtained.‘ In similar fashion, 
adiabatic compressibilities of a number of aque- 
ous ionic solutions,® and of pure water over a 
great range of temperature,*® have been precisely 
measured. 





*Submitted in partial fulfillment of the requirements 
for the degree of Doctor of Philosophy in the Graduate 
College of the State University of Iowa, Iowa City, Iowa. 

** Now of Sylvania Electric Products Inc., Flushing, 
Long Island, New York. 

'E. H. Amagat, Comptes rendus 104, 1159 (1887). 
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?R. E. Gibson, J. Am. Chem. Soc. 59, 1521 (1937); 
A. F. Scott, V. M. Obenhaus, and R. W. Wilson, J. Phys. 
Chem. 38, 931 (1934); D. Tyrer, J. Chem. Soc. 105, 2534 
(1914). 

4C. Bachem and E. Hiedemann, Zeits. f. Physik 94, 68 
(1935); E. B. Freyer, J. C. Hubbard, and D. H. Andrews, 
J. Am. Chem. Soc. 51, 759 (1929); S. Parthasarathy, 
Proc. Ind. Acad. Sci. 2A, 497 (1935). 

5C. Bachem, Zeits. f. Physik 101, 541 (1936); E. B. 
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APPARATUS AND GENERAL METHOD 


The method adopted for measurement involves 
the generation of a supersonic wave by a driven 
quartz crystal, and wave-length measurement by 
the optical method of secondary interferences.’ 
The commercial x-cut quartz crystal plate, about 
2.5 cm square and 2 mm thick, was driven 
indirectly by a radiofrequency power oscillator 
controlled by a quartz crystal. A single RK-20A 
tube was used in a conventional pentode circuit. 
The power supply was well filtered and regulated, 
and was more than adequate. A condenser-tuned 
resonant excitation control gave certain decided 
advantages: (a) control was smoothly variable 
over a wide range while leaving the power 
oscillator in operation under optimum condi- 
tions; (b) link coupling effectively insulated the 
control circuit and the crystal from the high 
supply voltages required by the power oscillator 
and permitted a convenient physical separation ; 
(c) use of a current indicator permitted repro- 
ducible excitation values. 

The driven quartz crystal was shellacked to a 
brass plate, which was mounted vertically within 
a brass box about 5.2 cm wide, 3.7 cm thick, and 
6.5 cm high, to be itself submerged in the liquid 
being measured. The grounded supporting plate 
and a tinfoil sheet shellacked to the front face of 
the crystal supplied its electrodes. The emitting 
front surface of the crystal was made closely 
parallel to the plane front face of the crystal 


7C. Bachem, Zeits. f. Physik 101, 541 (1936); C. 
Bachem and E. Hiedemann, Zeits. f. Physik 94, 68 (1935); 
N. Seifen, Zeits. f. Physik 108, 681 (1937-1938). 
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box, and lay about 6 mm directly behind the 
circular 2.4-cm aperture in the box face. The 
aperture was sealed with a thin mica diaphragm, 
which retained the xylol with which the crystal 
box was filled. 


Actual measurement was made in a wave 
trough about 40 cm X 7 cm X7 cm with the 
crystal box at one end and a plane brass reflector 
at the other. An interlocking brass framework 
around the open top of the trough provided 
pivot supports for the reflector and sliding ways 
for the support of the crystal box above the 
parallel plate glass sidewalls. Reflector orienta- 
tion and crystal box position and orientation 
were controlled by propelling screws opposed by 
springs under tension. With the wave trough 
nearly filled with the liquid, with the driven 
crystal properly excited, and with correct orien- 
tation and separation of reflector and crystal 
box, standing waves were set up in the liquid, 
having a direction of propagation along the 
trough length. 

Because standing density variations exist, the 
supersonic wave can be made visible by sending 
through it a transverse parallel beam of light, 
which is diffracted in very well-known fashion.® 
After emerging from such a diffracting medium, 
the light beam develops intensity variations® 
(“secondary interference’), which can be seen 
directly upon a screen placed in its path. The 
time average of the pattern as observed by the 
eye or recorded upon a photographic plate is 
that of the “sound grating’ producing the 
diffraction : narrow vertical bright lines separated 
by the half wave-length (A/2) of the supersonic 
wave in the liquid, or about ? mm in this 
research. 

From a 110-volt, 15-amp., d.c. carbon arc, a 
condensing lens, and an adjustable vertical slit, 
a divergent light beam traversed the wave 
trough about 50 cm distant. Some 50 cm beyond 
the trough a small, shielded lens of focal length 
approximately 4.5 cm cast a magnified image of 
the light pattern on a screen about 20 cm away. 
Separation and orientation of crystal box and 
reflector were carefully adjusted to give a simple 


§L. Bergmann, Ultrasonics (John Wiley & Sons, New 
York, 1938), English trans., p. 137; J. C. Hubbard, Am. 
J. Phys. 8, 207 (1940); C. V. Raman and N. S. N. Nath, 
Proc. Ind. Acad. Sci. 2A, 406 and 413 (1935). 

*N.S. N. Nath, Proc. Ind. Acad. Sci. 4A, 262 (1936). 
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sharp pattern of vertical bright lines. Figure 1 
gives a diagrammatic sketch of the optical 
system. 

The wave trough was mounted on the sliding 
bearing-block of a Gaertner precision micrometer 
screw of half-millimeter pitch, which permitted 
a maximum excursion of about 6 cm. Displace- 
ments were read to 0.005 mm. When the trough 
was displaced in the direction of its length, the 
entire pattern moved across the screen with a 
spacial period of 4/2. It was found convenient 
to measure seventy half wave-lengths of the 
supersonic wave for each liquid. Such a method 
of measurement did not require parallel light, 
and permitted additional magnification of the 
light pattern. 

The operating frequency was first measured 
by a General Radio Type 224 precision wave 
meter of the resonance type, which indicated 
approximately 1016 kc/sec. Then with a radio 
receiver tuned to KDKA, Westinghouse Station, 
Pittsburgh, Pennsylvania, an audible beat note 
was generated between the signal from the power 
oscillator, operating under load, and the carrier 
of KDKA. The frequency of the heterodyne 
represented the difference in frequency of the two 
radiofrequency waves and remained constant 
regardless of modulation of the commercial 
carrier.” Heterodyne frequency was measured 
as 4,760 c/sec. with a possible error of about 
10 c/sec. Since KDKA maintains a carrier fre- 
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quency within 2 c/sec. of its assigned 1020 
kc/sec., the precise frequency used in this 
research was 1015.24 kc/sec., with a possible 
error of about 10 c/sec. at 25°C. Room tempera- 
ture variations in the course of measurement 
were not important in introducing frequency 
change. 

From measured values of frequency and of 
distance D=70(A/2) cm for the liquid used, the 


1° W. Bragg, J. Acous. Soc. Am. 11, 5 (1939). 
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TABLE Ia. Summary of ccc measurements 
at 25°C. 














Am 
V(m/s) 10°98 =10"48 X10" am 
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Substance Mol.wt. M.F 
NH.Cl 53.497 .01999 1543.31 41.389 3.329 1.665 1.18 
LiNO; 68.948 .005471 1509.17 43.512 1.206 2.204 1.6 
LiCl 42.397 .01689 1551.14 40.767 3.951 2339 1.38 
NH.C2H:02 77.064 .01040 1534.46 42.214 2.504 2.408 1.00 
KCN 65.104 .009193 1525.76 42.340 2.378 2.587 73 
KCl 74.553 000056 1526.05 41.993 2.725 2.737 2.25 
NaCl 58.454 .01998 1564.77 39.230 5.488 2.747 1.48 
KCIO; 122.553 .003075 1502.27 43.861 0.857 2.787 8.0 
CoSO, 155.004 .002395 1509.95 43.861 0.857 3.578 7.3 
NaOH 40.005 .02090 1610.32 36.839 7.879 3.770 2.0 
Ca(NOs)2 164.086 .003937 1509.81 42.875 1.843 4.681 2.2 
Mg(NOs)2 148.336 .004041 1514.30 42.702 2.016 4.989 4.6 
Cle 95.234 006012 1532.72 41.636 3.082 5.126 38 
Co(NOs)2 182.956 .002006 1502.85 43.688 1.030 5.135 36.00 
iCle 129.604 004117 1516.04 42.500 2.218 5387 14.0 
CoCle 129.854 .003289 1512.85 42.922 1.796 5.461 8.1 
(NHs)280, 132.144. 1523.88 42.469 2.249 5.617 5.3 
uSO. 159.634. 1512.56 42.344 2.374 6.016 9.4 
AICls 133.341 .004656° 1533.59 41.595 3.123 6.707 7.9 
MgSOx 120.384 .003107 151842 42.619 2.099 6.756 7.5 
K2Cr207 204.212 .001984 1504.30 43.368 1.350 6.804 16.0 
Na2CO; 105.994 .01006 1589.14 37.531 7.187 7.144 No inf. 
Al(NOs)3 212.994 .001950 1509.14 43.268 1450 7.436 6.2 
K2S0. 174.256 .008918 1526.43 41.798 2.920 7453 9.5 
Co(CoHs02)2 176.988 001754 1510.74 43.390 1.328 7.571 60.0 
K:CO; 138.192 .003857 1530.11 41.757 2.961 7.677 4.6 
Ba(OH): 171.386 008125 1516.77 42.290 2428 7.770 23.1 
H:0 18.016 1497.62 44.718 
Tap water 1497.62 44.703 
raBLE Ib. Summary of paney measurements 
at 2 
AD/At Ao/. AD/at Av/At 
Substance X1e Am «my /°C) Substance X10 (cm/sec./°C) 
NHAC 9.54 277 CoCle 8.77 254 
LiNOs 8.92 259 (NH4)280, 8.91 258 
LiCl 7.39 “214 CuSO, 8.00 232 
NH.C:H;0: 8.00 232 AICls 8.31 241 
KCN 8.00 232 MgSO, 9.08 263 
KCl 10.77 312 K 9.23 268 
NaCl 7.08 205 Na:CO; 8.46 245 
KCIO; 9.69 281 Al(NOs)s 9.54 277 
CoSO, 11.23 326 r 9.08 263 
NaOH 6.77 196 Co(C2H302)2 8 62 250 
Ca(NOs)2 10.40 302 K:CO; 10.31 299 
Mag(NOs)2 9.23 268 Ba(OH): 8.77 254 
MgCle 7.39 214 
Co(NOs)s 8.92 259 H:0 9.54 277 
NiCh 9.39 272 Tap water 9.54 277 











supersonic velocity is found from the relationship 
v=vdr to be v=0.029007D X10*® cm/sec. Meas- 
urements by many investigators" show no cer- 
tain dispersion of acoustic velocity in liquids for 
frequencies below about 10° c/sec. 

Liquid density p was measured directly with 
a 25-cc pyknometer and sensitive analytical 
balances. Since the adiabatic compressibility 
8=1/pv’, fractional error in 8 can be expressed : 


Ag Ap 2Av [2 nat ~) ; 
















B p ov p 


1A. K. Dutta, Physik. Zeits. 39, 186 (1938); H. S 
Kneser, Ann. d. Physik (Leipzig) 37, 277 (1938); 
Parthasarathy, Proc. ind. Acad. Sci. 4A. 17 (1936); B. v 
R. Rao and D. S. S. Ramiya, Phys. Rev. 60, 615 (1941); 
B. Spakovskij, Comptes rendus Acad. Sci. U.R.S.S. (N.S.) 
18, 169 (1938). 
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The quantity AX/A is the largest contributor to 
AB/B and has a value less than 8X10~. 


EXPERIMENTAL PROCEDURE. RESULTS 


Solutions of the selected twenty-seven electro- 
lytes were made with doubly-distilled water.” 
Each was at a single concentration ranging 
from a mole fraction of 0.02 for small molecular 
weights to 0.002 for the highest molecular 
weights used. Each solution, of approximately 
1.6-liters volume, was carefully protected in a 
glass bottle with paraffin-coated cap. Contami- 
nation of the solutions during the course of 
measurement was largely prevented by coating 
with paraffin all surfaces excepting the glass 
sidewalls of the wave trough and the mica 
window of the crystal box, and by avoiding 
personal contact. 

Temperature gradients in the liquid were 
effectively prevented by a motor-driven glass 
screw 2 cm in diameter and 8 cm long, operated 
at 4 or 5 r.p.s. in the wave trough outside the 
supersonic beam. However, no attempt was made 
to insulate the trough thermally or to preserve 
a constant temperature. Cooled at first to about 
20°C, the liquid was allowed to change smoothly 
in temperature, during which time a number of 
quick determinations of wave-length and of 
density were made. Linear experimental curves 
of D=70(A/2) against ¢ and of p against ¢ were 
then constructed. From these curves, v and 8 
were computed at 25°C. Measurements of 
temperature were made by an immersed ther- 
mometer. 

Tables Ia and Ib present a summary of 
measurements at 25°C and a number of quanti- 
ties obtained by further computation. In order 
to avoid error due to water of crystallization and 
moisture content of the solute crystals, P, the 
percent by weight of the solute, was determined 
for each solution by linear interpolation at a 
temperature for which both P and p are given in 
International Critical Tables, Vol. III. M.F. is 
the mole fraction of the solute. The symbol Ag 
signifies (Bwater —Asolution) at the temperature of 
measurement. The decrease in 8 per mole fraction 


2 The doubly-distilled water, prepared by W. W. 
Boardman of the State University of lowa Chemistry 
Department, had an approximate conductivity of 5x 10- 
mho/cm. 
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of the solute is given the symbol —A6/Am. The 
average increase in sound velocity per degree 
centigrade for each solution over the temperature 
range of 23°C to 26°C is given in Table Ib. 

With any method of measuring sound velocity, 
excessive sound amplitude may cause errors. A 
high temperature gradient may exist in front of 
the emitting crystal. The emitted wave may 
not be plane, and in interferometric methods of 
measurement, successive positions of maximum 
reaction, chosen on the basis of average effect of 
all components, are not necessarily separated by 
the half wave-length." It has been reported that 
the supersonic intensity required for velocity 
measurement with the supersonic interferometer 
for liquids is significantly less than the intensity 
required by optical methods ;“ but as can be seen 
from Table II, the sound intensity employed in 
this research has not been great enough to 
impair accuracy of measurement. 

Within the limits of experimental errors, sound 
velocity is measured to be the same in tap water 
as in the most carefully prepared distilled water 
(Table Ia). Additional tests showed no velocity 
difference when ordinary distilled water was 
used, whether or not it had been freed of dis- 
solved air by boiling. These results are in ap- 
parent disagreement with reports!® which indi- 
cate a strong dependence of velocity upon gas 
content. It would seem that the discrepancy 
arises from two sources. First, experimental 
conditions under which the gas appears as fine 
bubbles throughout the liquid have not been 
carefully distinguished from conditions under 
which no bubbles appear; the liquid is clear and 
homogeneous; and the gas is truly dissolved. 
Second, since agitation of the liquid by the 
acoustic wave aids any tendency of the gas to 
come out of solution and to form bubbles, 
methods of measurement are subject to error if 
the sound intensity employed is so great that 
bubbles of the once-dissolved gas are formed, or 
actual cavitation of the liquids is produced. 

In the present series of experiments the liquids 
. remained homogeneous and clear. There was no 


4 W. T. Richards, Rev. Mod. Phys. 11, 44-48 (1939). 
“W. Garten, Bull. Am. Phys. Soc. 16, 7, Abstract 32 
1941). 

%K. Dérsing, Ann. d. Physik 25, 228 (1908); Ch. 
Sérensen, Ann. d. Physik (Leipzig) 26, 121 (1936); Taylor 
and Sproule, Trans. Roy. Soc. Can. 3, 23, 91 (1929). 
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apparent formation of gas bubbles in the body 
of the liquid due to excessive sound intensity. 
Under such conditions it is not surprising that 
no noticeable change in sound velocity was 
observed, since even water saturated with air at 
25°C contains only two Nz molecules and one O, 
molecule to approximately 2105 molecules of 
water. Tap water from soft water sources may 
likewise have a low ionic concentration and an 
adiabatic compressibility which differs only 
slightly from the value for pure water. 


DISCUSSION, INTERPRETATION, CORRELATION 
OF RESULTS 


Precision compressibility measurements by 
Randall'® have shown that with increasing 
temperature, the adiabatic compressibility for 
pure water decreases to a minimum about 64°C, 
then rises again. Not long after the early work 
of Amagat, Tamman" ennunciated a theory 
attributing the decreased compressibility of 
aqueous solutions to increased internal pressure. 
On the same basis, since increased external 
pressure was known to lower the temperature of 
minimum compressibility for water, he antici- 
pated the same behavior for aqueous solutions 
due to increased internal pressure and supported 
the theory for the case of ethyl alcohol. Tam- 
man’s principle may be stated: ‘‘A given amount 


TABLE II. Comparison of measured values of v and £8 for 
pure water. 








Freq. 
(ke/sec.) 
1015.24 
1015.24 


Observer 
Corey 
Corey 
Seifen 
Schaaffs 
Bachem 


Hiedemann 
Seifen and 


Schreuer 
Parthasarathy 


10'°8 
44.718 
44.872 


. v (m/sec.) 
1497.62 
1494.87 
1497.15 
1467 
149741 
1497.18 





interference 
Secondary 
interference 


1494 J Deb 
grating 
1465 Deb: 
grating 
1497.4 Interferometer 
1 


483.1 Interferometer 
1495 Interferometer 


Bergmann 


Randall 
oy 

reyer, 
Hubbard, 
Andrews 
Hubbard 1498.1 Interferometer 
and Loomis 


Tyrer Piezometer 








16 See reference 6. 
17G. Tamman, Zeits. f. physik. Chemie 11, 676 (1893). 
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Fic. 2. Correlation between —10A8/Am and —dah/am 
for the following numbered electrolytes: 1. NH,Cl, 2. LiNOs, 
3. LiCl, 4. NH«C2HsO2, 5. KCN, 6. KCI, 7. NaCl, 8. KCIOs, 
9. CoSQ,, 10. NaOH, 11. Ca(NOs)2, 12. Mg(NOs)2, 13. 
MgCl, 14. Co(NOs)2, 15. NiCle, 16. CoCle, 17. (NH«4)2SO,, 
18. CuSQ,, 19. AICI, 20. MgSO, 21. KeCr2O;, 22. Al(NOs)3, 
23. K:SO, 24. Co(C2H;02)s, 25. K:COs,, 26. Ba(OH)>. 


of water in a solution behaves like the same 
amount of pure water subjected to some greater 
external pressure.’”’ A special series of measure- 
ments to 40°C on pure water and on a 5.6 percent 
solution of Na2CO; gave results which appeared 
to be in accord with Tamman’s principle, but 
showed definitely that even for maximum solute 
concentrations used, the compressibility mini- 
mum for the solution lay above even 50°C. 
Bachem!*® and Freyer” have shown that the 
compressibility of an aqueous ionic solution 
decreases only slightly less than linearly with 
solute concentration over a great concentration 
range. Therefore, since concentrations used in 
the present research are 0.01 mole fraction or 
less, with a few exceptions of approximately 0.02 
mole fraction, values of —10"A8/Am as they 
appear in Table Ia give a valid comparison of 
the behavior of different electrolytes. Since 
values of —10"A8/Am for the solutions are more 
or less characteristic of the solutes to which they 
pertain, a correlation is expected between these 
values and other mole fraction indications of 
change, such as a change in water structure as 
‘ given by x-ray experiments. X-ray diffraction 
curves exhibiting maxima and minima, such as 
those for water and other liquids, may be 
produced by the existence of a few fairly definite 
interatomic distances.” From an experimental 


18 See reference 5. 
19 E, B. Freyer, J. Am. Chem. Soc. 53, 1313 (1931). 


20 B. E. Warren, J. App. Phys. 8, 645 (1937). 
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x-ray diffraction curve for a liquid, a unique 
radial distribution function may be directly 
determined, which yields information regarding 
number and arrangement of neighboring mole- 
cules, or the structure of the liquid. Changes in 
structure under different experimental conditions 
may thus be examined. But as might be expected, 
the diffraction curve itself shows an associated 
change. In the case of the x-ray diffraction curve 
for water, the second peak is particularly sensi- 
tive to changes in the radial distribution func- 
tion, and by its change in height may serve as 
an indicator of the effects produced.” Introduc- 
tion of ions, like an increase in temperature, has 
in general the effect of breaking down the 
structure and causing closer packing; and the 
second diffraction peak shows a correlated de- 
crease in height as the coordination number 
increases. The quantity —dh/dm, proportional 
to the change in height of the peak per mole 
fraction of the solute, is an indication of the 
structural change characteristic of the solute.” 

Figure 2 shows graphically the correlation 
between —0h/dm and —10”"A6/Am for twenty- 
six of the twenty-seven electrolytes considered 
in this research. Most of the observations lie 
within a rather small angle. The grouping of the 
univalent electrolytes, especially of LiNOs, LiCl, 
NH,C:H;0O2, KCN, KCl, and NaCl, indicates 
that the observed effects of an ion depend much 
more strongly upon its charge than upon its 
radius. 

In general, a greater effect results from electro- 
lytes of higher valence types. It is clear, however, 
that the effects of individual ions are not easily 
distinguishable. The behavior of an ion is much 
influenced by the character of the ion with 
which it is paired. For this reason, although one 
may arrange a sequence which is valid for a 
number of cases in a variety of experiments, one 
cannot set down unequivocally a series of cations 
or of anions in order of their effects. 

A comparison of results shows satisfactory 
agreement with: conclusions of Passynski,™ 
whose work was based upon a theory of Debye; 


#1 J. Morgan and B. E. Warren, J. Chem. Phys. 6, 666 


(1938); G. W. Stewart, Phys. Rev. 37, 9 (1931); G. W. 
Stewart, J. Chem. Phys. 7, 869 (1939). 

2G. W. Stewart, J. Chem. Phys. 7, 869 (1939). 

%3 A. Passynski, Acta Physicochimica (8) 4, 385 (1938). 

*“P. Debye, Festschrift H. Zangger (Rasch and Co., 
1935), Part II, 877. 
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the effects expected according to the theory of 
‘permanent hydration,’ due to Bernal and 
Fowler ;?5 the influence of ionic radius as given 
by Pauling;?* and the evidence of change in 
water structure found by Frank and Robinson?’ 
in computed partial molal entropies. Because to 
a first approximation an ionic solution of low 
concentration may be looked upon as altered 
water, one expects a correlation between prop- 
erties of the solution and properties of the water 
within the solution, i.e., the solvent water. 
From the same point of view it is of interest to 
compare at 25°C the values of solution compressi- 
bility, 8, with corresponding values of V;, the 
partial molal volume of the solvent water. By 
use of the developments of Geffcken®* and 
Wirth,” an expression has been derived for V; 
which is applicable in all cases where ®, the 
apparent molal ionic volume of the solute, varies 
linearly with some power of the solute concen- 
tration. For sixteen of the twenty-seven electro- 
lytes used in this research, @ does not vary 
linearly with c! (cf. Stewart®) as the simple 
Debye-Hiickel theory predicts. The derivation 
of the expression follows immediately. 


Explanation of symbols used: 


V: volume of the solution in milliliters; 

n: moles of designated substance in the solution; 

c: moles of designated substance per liter of solution; 

V=2aV/dan: partial molal volume of designated substance; 

#(V)=¢: apparent molal ionic volume of designated 
substance; 

x: exponent of ¢ which makes the plot of ® against c* a 
straight line; 

k: slope of linear plot of @ against c*; 

Subscript 1 designates the solvent; 

Subscript 2 designates the solute; 

Subscript or superscript 0 designates the value of a 
quantity when the solute concentration is zero. 


If temperature and pressure remain constant, 
a small change in volume may be expressed 
generally as 


dV=(aV/am;)dn;+(aV/dns)dnz. —(1) 


J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 
515 (1933). 

1. Pauling, The Nature of the Chemical Bond (Cornell 
University Press, 1940), pp. 47, 348. 

27H. S. Frank and A. L. Robinson, J. Chem. Phys. 8, 
933 (1940). 

28 W. Geffcken, Zeits. f. physik. Chemie 167A, 240 (1933). 

29H. E. Wirth, J. Am. Chem. Soc. 62, 1128 (1940). 
% G. W. Stewart, J. Chem. Phys. 7, 381 (1939). 
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If the mole ratio also remains constant: 


dV/V=dn,/n,=dn2/n2=a. (2) 


Thus, upon substitution: 
V =n,(0V/dn;)+n2(0V/dnz2). (3) 


Since V,;=(0V/dm), etc., and co=1000n2/V, 
etc., 
V—n2V2 1000-2, 


i= 





(3a) 
nN) Ci 


By definition, 
V—Vo 





(V2) =%,= (4) 


Ne 
From (3), Vo=m:V1°; therefore 
1000 —c,V;° 


d 1000 — doc, 
@, =——__—_——__,. an G4 Rowe, 
Ce V,° 


Therefore from (3a) 
Vi =(V 191000 —c2V2)/(1000—c22). (3b) 
From (4), V = Vo-+m2%2; therefore, 
(8V/dn2) = V2 =22-+02(db2/dn2) 
= 2-+12(dB2/dC2)(9c2/dn2). (5) 
With the exponent x chosen empirically, the 


plot of 2 against c* may be made sensibly linear 
e — 





TABLE III. Summary of measurements and computations 
on 6 and V;. The temperature in each case is 25°C. 











Substance Exp. of V» k ce: Vv, 10348 
NH.Cl 0.5 36.30 1.225 1.0827 18.05642 41.389 
LiNOs 0.8 28.99 605 0.3018 18.06839 43.512 
LiCl 0.5 17.02 1.46 0.9342 18.05727 40.767 
NH.C3:H30: 1.0 60.72 .74 0.5616 18.06503 42.214 
KCN 1.0 30.5 .22 0.5055 18.06838 42.340 
KCl 0.5 26.36 241 0.5480 18.06044 41.993 
NaCl 0.5 16.28 2.22 1.1048 18.04563 39.230 
KC1O;3 1.0 45.78 —4.87 0.1694 18.07196 43.861 
CoSO. 10 —-19 15.0 0.1329 18.06463 43.861 
NaOH 10 -—3.75 1.47 1.1840 18.03226 36.839 
Ca(NOs): 0.8 42.2 5.08 0.2167 18.06469 42.875 
Mec 1.5 39.57 3.29 0.2226 18.06731 42.702 

I: 0.5 15.8 5.88 0.3326 18.05916 41.636 
Co(NOs): 0.5 28.7 13.13 0.1108 18.06502 43.688 

iCls 0.5 16.42 —6.58 0.2281 18.07590 42.500 
CoCle 0.5 13.9 2.92 0.1821 18.06734 42.922 
(NH4)3SO4 08 53.35 6.61 0.2198 18.06307 42.469 

uSO« 05 —49 9.72 0.2193 18.06039 42.344 
AICls 0.5 12.7 9.76 0.2565 18.05790 41.595 
MgSO. 0.7 —1 8.29 0.1725 18.06413 42.619 
KeCrsOr 1.0 92.45 11.52 0.1089 18.06692 43.368 
NazCO; 0.7 —3.9 9.78 0.5618 18.02279 37.531 
AI(NOs)s 0.5 47.30 4.93 0.1076 18.06783 43.268 
K2S04 1.0 27.4 6.09 0.2723 18.06118 41.798 
Co(C2HsO2): 1.0 34.0 101.1 0.09681 18.05222 43.390 
KsCOs 0.8 14.9 8.55 0.2135 18.06170 41.757 
Ba(OH): 1.0 2.19 —37.73 0.1737 18.08996 42.290 
H:0 18.06941 44.718 
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over a great concentration range. With x so 
_ chosen, the equation ©, = 6,°+c,? is linear, with 
slope k and intercept #,°. 
Since experimentally 2. = 9+ ke2", 


09,/dc2=kxc2*—"', 
c2=1000n,2/ V, 
V—n2(0V/dn2) 


OC2/ON2 = 1000 
V2 





From (5) with substitutions: 


OV B® + R(x +1) C27 
—_—_= V2 = 


ang 1+ (kxcs*+1/1000) 





As one can easily verify from columns 3, 5, 
and 6 of Table III, the second term in the 
denominator of (6) is negligible, especially in 
comparison with errors in empirical determina- 
tion of #,°, k, and x. Therefore with this approxi- 
mation it results that 


V2=8,9+k(x+1)c2". (7) 


That 4,°=V,° is of course apparent from defini- 
tions of $2 and V2. 
Substituting (7) into (3b) there results: 


1000 — co2°— k(x+ 1)c,7*! 


Vi= V;° 
1000 — coP2° — keg? *! 


(8) 





When %, varies linearly with c (i.e., x=4), 
(8) reduces to Wirth’s formula for V}. 

Equation (4) which defines @(V2) is readily 
transformed to an expression suitable for compu- 
tation. Introducing the density of the solution p, 
the molecular weight of the solute M2, and the 
density of the solvent p, it follows that: 


“ Mz 1000 (p—p1) 
(V2) =%,=—— ; 
P1 Pi C2 





C2=10Pp/M, (P=percent by weight of solute). 

Table III was prepared on the basis of the 
preceding development. From information given 
in the International. Critical Tables, Vol. III, 
supplemented in a few cases by density-tempera- 
ture data from this research, values of #2 and 
corresponding: values of c, were computed at 
25°C. For each electrolyte a sufficient number of 
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Fic. 3. Correlation between adiabatic compressibility 
of solutions and partial molal volume of solvent water at 
25°C for the following numbered electrolytes: 1. NH,Cl, 
2. LiNOs, 3. LiCl, 4. NH«C2H3O2, 5. KCN, 6. KCI, 7. NaCl, 
8. KCIO;, 9. CoSOQ,, 10. NaOH, 11. Ca(NOs)2, 12. Mg- 
(NOs)2, 13. MgCle, 14. Co(NOs)2, 15. NiCls, 16. CoCls, 
17. (NH4)2SO,, 18. CuSO, 19. AlCl, 20. MgSO, 21. 
Ke2CreO7, 22. Al(NOs)s, 23. K2S0O,, 24. Co(C2H;02)s, 25. 
K;CO;, 26. Ba(OH)>. 


plots of @, against powers of c; were made to 
permit the selection of the linear graph. Slope, 
intercept, and exponent as determined in this 
manner appear directly in Table III, and as well 
were used together with proper values of Cc» in 
finding V; from (8). It is true that the values 
of &,° appearing in Table III have been obtained 
by extrapolation and are not necessarily the 
values which ®, would assume at zero concen- 
tration. However, a comparison of (3b) and (8) 
reveals that ,° and V,° appear in (8) only as 
part of &, and V2, respectively; and that conse- 
quently, linearity of the plot of 2 against c* in 
the concentration range containing the measure- 
ment is the sole requirement for the applicability 
of (8). The interpolation factor of (8) has 
virtually the same value whether the cubic 
centimeter or the milliliter is used as the volume 
unit. Thus the volume unit in terms of which V; 
is expressed in Table III depends upon the fact 
that V,° is expressed in cubic centimeters. 

The circles of Fig. 3 represent experimental 
information on the electrolytes of this research, 
each at a single concentration (Tables I, III). 
The dots defining the two curves represent values 
of B at solute concentrations of 0, 0.1, 0.25, 0.35, 
0.5, 0.65, 0.8, and 1.0 mole per liter, respectively, 
found by graphical interpolation from Bachem’s 
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work,*! and corresponding values of V; computed 
therefrom by the author in the manner already 
described. In most cases, irrespective of the 
type of solute, there is a surprising correlation 
between corresponding values of 8 and V. 

The significance of the correlation is more 
readily appreciated through physical interpreta- 
tion. In nearly all cases the addition of an 
electrolyte produces two intuitively compatible 
results: the partial molal volume of the solvent 
decreases, and the compressibility of the solution 
becomes less. These are the results one would 
anticipate if the water structure were more 
highly coordinated by the addition of ions. On 
the other hand, if compressibility values in- 
creased in general with decreasing partial molal 
volume of the solvent, one could not confidently 
maintain that a more highly coordinated water 
structure is responsible for the observed change. 

That the two full curves from Bachem diverge 
slightly is not at all strange when one considers 
the effects which such factors as ionic radius, 
configuration, associated ion pairs, etc., might 
have upon both 8 and V;. More significant is the 
surprisingly close correlation which is shown 
between these two curves, and by the general 
clustering of the circles for nearly all of the 
electrolytes investigated. One is led to the 
conclusion that in most cases, solvent structure 
change exercises an unmistakable influence upon 
both 6 and V. 

According to the International Critical Tables 
from which values of 2, are computed, KCIQOs, 
NiCl:, and Ba(OH): are peculiar in that their 
apparent molal ionic volumes decrease with 
concentration. The slope of the curve for NiCl. 
changes, and even reverses in sign as the temper- 
ature rises from 17°C to 25°C. These cases, 
apparently anomalous, are represented on Fig. 3 
by circles 8, 15, and 26. Curves belonging to 


3! See reference 5. 
® Circle 26 is outside graphed region. 
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these electrolytes would thus show an increasing 
8 with decreasing V;, in contrast with all other 
electrolytes studied. Recent measurements by 
Redlich and Bigeleisen,* however, have shown 
the Critical Tables to be in error with regard to 
the anomalous character of KCIQs. It is possible, 
therefore, that the Tables are erroneous in the 
information which predicts anomalous behavior 
of NiCl, and Ba(OH)s. 

A perfect correlation between 6 and JV; is 
hardly to be expected since there are a number 
of contributing factors in each case and since 
computation of V, involves an empirical evalua- 
tion of constants. The reason is not obviously 
indicated for the peculiar position of Co(C2H3O2)s, 
which seems to have more effect upon solvent 
structure than upon solution compressibility 
(cf. Fig. 2); or for the possibly anomalous 
behavior of NiCl, and Ba(OH)s, for which the 
curves of ®. against c* have negative slopes at 
25°C. In such unusual cases there are clearly 
other agencies besides the compacted solvent 
structure which account for the observed effects. 

From a positive standpoint, however, the 
results are definite and heartening. For all 
electrolytes the mole fractional change in com- 
pressibility is satisfactorily correlated with x-ray 
and other indications of liquid structure change 
in the water; and the compressibilities are them- 
selves correlated in a striking manner with 
corresponding values of the partial molal volume 
of the solvent water. All such facts strongly 
support the conclusion that the water structure 
becomes more highly coordinated and compacted 
with the introduction of ions. No _ well-sub- 
stantiated evidence presented in this paper 
maintains a contrary view. Further research has 
been instituted on some of the problems which 
have arisen during the prosecution of this 
research. 


3% Redlich and Bigeleisen, J. Am. Chem. Soc. 64, 758 
(1942). 
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In this paper, the Schrédinger equation is solved analytically for the Coulomb potential 
with horizontal cut-off. For this purpose, a detailed mathematical study of the confluent 
hypergeometric functions is carried out and their similarity to Bessel and Neumann functions 
is put on a quantitative basis. The resulting formulas for the energy levels are valid for all 
cases except when the energy is lying very little above the cut-off. 
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HE Coulomb potential with a cut-off (see 

Fig. 1) has been proposed as an approxi- 
mate potential energy for certain quantum 
mechanical problems in the theory of solids.! 
It has actually been used for such a purpose by 
Tibbs, who solved the problem in one particular 
case by numerical integration methods.’ But, 
in fact, the ‘problem is capable of an analytic, 
though not rigorous, solution which is to be 
derived in the following. 

The problem is to solve the equation 


ad°F a 2a I(l+1) 
—+(-<+=-——)r-0 
dr? 


k? ro r? 
















for r<ro and 









d?F a? 2a Ki+1) 
—,+(-— — )F=0 (2) 


dr? 


« rT r2 





for r>ro, and to join them with continuous 
derivative at r=ro. The usual conditions are to 
be added at the origin and infinity. The energy 
parameter is written here in the form —a?/k’, 
where a is a given constant and k is our unknown 
to be determined. 










1. BESSEL-LIKE AND NEUMANN-LIKE FUNCTIONS 





The solutions of Eq. (2) are called confluent 
hypergeometric functions and are discussed in 
Whittaker and Watson.’ The two solutions may 
be either distinguished by their behavior at the 
origin or their behavior at infinity. In the nota- 
tion of WW, we have one solution, M,, 14.;(2ar/k), 





1G. H. Wannier, Phys. Rev. 52, 191 (1937). 

2S. R. Tibbs, Trans. Faraday Soc. 35, 1471 (1939). 

3E. T. Whittaker and G. N. Watson, A Course of 
Modern Analysis (Cambridge University Press), Chapter 
XVI. In the ifowlor quoted as WW. 
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which is finite at the origin while the other, 
M,,-1-;(2ar/k), has a pole of order /. At plus 
infinity, one solution, W,,14;(2er/k), is finite 
while the other tends to infinity. It is well 
known that whenever the expression k—/ is a 
positive integer, the solution which is finite at 
the origin also happens to be finite at +, 
Hence k =/+1,/+2,/+3, --- are the eigenvalues 
of k in the case of the Coulomb potential without 
cut-off. If there is a cut-off at ro, then Wg, 14,(2) 
is the correct wave function beyond fo, but the 
values of & will no longer be integer. 

The similarity between the functions M and 
the Bessel functions is well known.‘ In order to 
bring it out we introduce a new variable x 
which is defined thus: 


x= 2(2ar)'=2(kz)! (3) 
and instead of M we introduce J,*(x), 


zip 


p! 





J,*(x)= My, ip(2). (4) 


Ee Eee 





f 


| 

} 
| 
| 


Fic. 1. Coulomb potential with cut-off. Solving the 
Schrédinger equation for such a potential is the problem 


of this paper. 


4H. A. Bethe, Handbuch der Physik, second edition, 
Vol. XXIV, 1, p. 287. 
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It can be verified by substitution into Eq. (2) 


that 
lim J,*(x) = J p(x), 


i.e., the zero energy Schroedinger eigenfunction 
is a Bessel function. In addition, by comparing 
power series, one finds that the first two terms 
of the series agree for all values of k. A better 
standard of comparison is gained by a study of 
their integral representations. The well-known 
Bessel integral may take the following form 


Ja)=— ff exp [=(«—-) Jorrru (5) 


The path comes from — ©, circles around the 
origin in a positive direction and returns to — « 
(Fig. 2). A similar integral exists for J,*(x). 
Whittaker and Watson introduce an integral 
representation for W;,, only.’ The same integral 
will also be valid for M;,» if the path is suitably 
altered. Transforming their integral with (3) 
and (4) and introducing a new variable of 
integration u through 


t= —}(xu—z2) = —4x[u+(x/4k)], 


we get 


1 x \~-A@tD 
IHay=— fen(ut—) 
2a1 4k 


x k—4(p+1) 
x («-—) du. (6a) 
4k 


The path in WW must be altered to encircle 
both singular points of the integrand in a 
counterclockwise direction and then to return to 
—« (Fig. 3). This statement is verified as 


j 1) 
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Fic. 2. Path of integration to be followed in the definition 


of the Bessel function by the Bessel integral [Eq. (5) ]. 





~ § WW, section 16-12. 
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Fic. 3. Path of integration to be followed in defining the 
Bessel-like function by a contour integral [Eq. (6). 


follows: If we substitute ¢ back again into (6) 
we get 


T= -ep ( )(as)°(1/2ni) 


x? 
8k 
x f e~*(—2)-?-1[1 + (x?/4et) PA@tds. 


Since the path encircles both singularities it can 
be deformed to make |t| >x*/4k. Hence the 
binomial can be developed into a power series, 
and the series can be integrated term by term. 
The resulting series, by (3) and (4), agrees with 
the series given for M in WW. 

To bring out the similarity between (5) and 
(6) we write the latter in the form 


J,*(x) = (1/2nt) 


1 
x f exp | 3x(»——) furrso(a/Aeupa, (6b) 
u 


with 





1—s}* 
¢(s) — (1 = stor em 
1+s 


=1+3(p+1)s*—fks*+---. 


Since |u| may be made larger than 1 in the 
integrand, we can conclude that 


x? 
J,(x) =Jp(x)] 140 8 
(x) (| + (—)] (8) 


(7) 





It is possible, incidentally, to derive a de- 
velopment of the Bessel-like function J,*(x) in a 
series of Bessel functions, following up the 
reasoning just indicated. It can be stated as 
follows: If we write for the series (7) 


$(s) =Lbs, (9a) 
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then 


@ 
J,*(x) = 2 b,(x/4k)’ J p+»(x). (9b) 
The convergence is quite rapid, but there are 
other, more physical ways of constructing im- 
proved expressions for J,*(x) which are outlined 
in Sections 3 and 4 below. 

In working with Bessel-like functions of integer 
index p we are faced with a difficulty which is 
well known from the theory of Bessel functions. 
Usually, J,* and J_,* form a complete solution 
of our differential Eq. (2), but whenever p 
becomes integer the two solutions become 
identical, apart from a constant factor. We must 
therefore define an auxiliary function to use 
beside J,*, and it is natural to look for a function 
which may be compared with the Neumann 
function N,(x) (called Y,(x) in some books). 
In the latter case, the definition runs 


N,(x) sin pr=J,(x) cos pr—J_,(x). (10) 


The right-hand side of this equation is con- 
structed so as to vanish identically for integer p. 
The same principle must be applied in con- 
structing a Neumann-like function N,* from the 
two Bessel-like functions J,* and J_,*. The 
relationship between the latter two may be 
obtained from their power series. If we add to 
this relationship two requirements which we wish 
to impose, namely, 








N,*(x) 
lim ———-=1 for all & 
7-9 N(x) 
and N,*(2) 
x 
lim =1 for all x, 
ke N,(x) 


then the new function is completely determined : 
N,*(x) sin pr 
_ (k+4p—9)! 
- (k-40-9) 





Jp*(x) cos pr—J_,*(x). (11) 


We can expect the new function to equal N,(x) 
under the same condition (8) as in the previous 
case. ® 

* It appears at first sight that the factor in front of Jp* 


alters the situation because it is possible to satisfy the 
requirement (8) 


x<K4k 


and yet to have that factor markedly different from 1 by 


The Bessel-like and the Neumann-like function 
form a complete system of solutions for (2), and 
therefore, our wave function W,,m(z) is expres- 
sible in terms of them. These relationships are 
well known.’ We may write the relation relevant 
to our case 


2 W,. m(z) sin 2mxr 
= (k-+-m—})!k-"Som (x) sin (k-+-m—})x 
—(k—m—4)!k"J-om(x) sin (k—m—4)e. (12) 
This relation reduces to 0=0 for integer or half 


. ° bos k ‘ 
integer m. But if we eliminate J_2, with the help 
of (11) we get a non-trivial relation for all m 


2? Wy, m(2) 
k 
= (k+m—})!k-"Jom(x) cos (R—m—})xr 
k ‘ 
—(k—m—}4)!k™Nom(x) sin (kR-—m—})xr. (13) 
Through its trigonometric factors, the formula 
shows up clearly the oscillating character of W 
as a function of k. Every time k—m++}3 is a 
positive integer the term in N will drop out. 
The oscillatory character does not continue, 
however, into the region where k <m because the 
vanishing of the sine is cancelled by the pole in 
the factorial. 


2. A SIMPLE APPROXIMATE SOLUTION 
OF THE PROBLEM 


Our problem consists in fitting the solution of 
(1), namely, 


x? Xo" i 
Fr=sJin(—(1- ) ). (14) 
2x0 16k? 


where we have defined consistently with (3) 


Xo = 2(2aro)! (15) 





and the solution of (2), which may be written 
according to (13), 


Fe=xJarai(x) cos kr 





(e-1-1) ee" ——* 
(kid! xNoi4i(x) sin kw. (16) 
having 
p~2k. 
But these two conditions imply also 
x<p 


and under these conditions, the lowest negative power in 
J_»* outweighs all other contributions to the function N. 
7™WW, section 16-41. 
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It is possible to get a quick and simple answer 
by observing the following. We have a condition 
under which the Bessel-like and Neumann-like 
function can be replaced by their namesakes, 
namely, 


x<Ak. (17) 


Now it is very fortunate that for physical 
reasons a somewhat weaker condition is always 
satisfied at the cut-off x =x: 


x <A4k. (18) 


We obtain the condition by demanding that the 
total energy be higher than the potential energy 
at the cut-off. Condition (17) will, in fact, apply 
to all but some of the lowest levels although the 
latter are no doubt the most important ones. 
Keeping this restriction in mind we shall con- 
struct an approximate solution valid for most 
states by setting equal at x =x» the logarithmic 
derivatives of the following functions 


Fy =xJ144(x?/2x0) (19) 
and 


F.=xJe141(x) cos kr 








(kR—1—1) !k?*! 
*Noiii(x) sin kw. (20) 
(k+/1)! 
We find the following result 
(k—1—1) !k2!+! 
—tan kr 
(k+l)! 


_ J 01(x)Ti44(3%) — Jor4s(%)Sia($x) 
Noi(x)Jiga(4x) — Neavss(x)Jia(4x) 


(21) 





where we have replaced for brevity the cut-off 
coordinate x» by x. Equation (21) is explicit in 
that it contains on the right a function of the 
cut-off coordinate only and on the left a function 
of the unknown k, the “effective” principal 
quantum number. For s states the formula is 
particularly attractive: 


Jo(x) sin $x—J;(x) cos $x 
No(x) sin $x —N,(x) cos hy 





(22) 


tan kr= — 


To realize the power of this simple formula we 
shall recompute with its help the Tesults of 
Tibbs.? Tibbs calculated in a particular case the 
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energy values for the 1s, 2s, and 2p states. 
His cut-off is located at 


x = 2.62 
and he found for k: 


k=1.154+.002 in the Is state, 
k=2.143+.019 in the 2s state. 


Formula (22) gives 
k=n+0.158 for all s states. 


The agreement may be partly accidental because 
it is upset if a refined approximation is used for 
one of the wave functions only. However, it is 
not impossible that formula (22) is actually 
better than it appears from this paper and that 
the deduction of the author is not doing it full 
justice. In any case, formulas (21) and (22) are 
very useful for a rough location of the energy 
levels. 

What follows first from (21) and (22) is that 
while the energy correction is of course smaller 
for high energy levels, the correction to the 
principal quantum number is of the same order 
for all levels having the same angular momentum 
quantum number. It is actually the same (within 
the restrictions of this section) for all s states 
and tends to be the same for sufficiently high 
energy in the other cases. If we turn to comparing 
states of different /, then s states are more 
disturbed than p states, and so forth. One can 
appreciate that difference particularly well if one 
assumes x sufficiently small to permit power 
series development of the Bessel and, Neumann 
functions. We get in that case from (21) 


(MD) 


($x) t+ . 
(21-+1) 1(21+3)! 





where n is the next lower integer, i.e., the 
principal quantum number in the correct sense 
of the word. The limiting formula for large x is 
also very simple if the Bessel and Neumann 
functions can be replaced by their trigonometric 
equivalents and if the ratio of factorials can be 
set equal to 1. We get 


k=n+(x/2nr) —}$1—}. 
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0.1 


x 





ie) 


Fic. 4. Plotting energy against cut-off radius according 
to Eq. on The abscissa x goes as the square root of the 
radius [Eq. (3)] and the ordinate contains the difference 
of the principal quantum number » and a number k which 
must be substituted into the Balmer formula to give the 
energy levels. The graph applies only to s states under 
the restriction (17). 


The validity of this formula is somewhat re- 
stricted because x must be large compared to /, 
yet condition (17) must not be broken. A 
graphical evaluation of formula (22), plotting 
k—n against x for s states is given in Fig. 4. 
One sees that the curve oscillates about the 
asymptote which we discussed just above. 


3. A MORE PRECISE SOLUTION FOR s STATES 


It has been pointed out that the solution of 
Section 2 is correct if condition (17) holds while 
only condition (18) is fulfilled in many actual 
cases. In this section, we shall improve Eq. (22) 
by applying various corrections to the approxi- 
mation (20), so as to make it quantitatively 
correct everywhere except in the immediate 
neighborhood of x = 4k. 

When one studies the approximate wave 
functions (20) numerically and compares them 
with the Schrédinger eigenfunctions in the case 
of integer k, one finds of course complete dis- 
agreement beyond x=4k. Below that point the 
two functions show the same general structure, 
but vary somewhat in detail. It is possible to 
cancel out most of this discrepancy by the 
following two steps. 

In the first step, we carry out a change in the 


abscissa scale so as to bring the nodes of the 
true wave function and its approximation to 
coincidence; i.e., we write the wave function in 
the form 


F=x,Z,(xx), 


where Z;(x) stands for the linear combination 
(20) of Bessel and Neumann functions. x; is an 
unknown function of x differing from it very 
little as long as condition (17) holds. If we 
introduce a corresponding r, through (3) we 
get for F the equation 


d*F 
(5;)+ )+ (= *)F=0 
while as a function of r, it obeys Eq. (2), 
d°F a? 2a 
ra | 
These two equations give for r, as a function of r 
ee 2a 2a / dr; 1 dF d’r, 
Ma? sae 


It is approximately correct to leave off the 
second term on the right-hand side because 
(1/F)(dF/dr,) is an expression of alternating 
sign which will make little long range contribu- 
tion to (dr,/dr)*. The remainder of the equation 
reads then 

dr; 1 a i 

aft 

dr , 


or passing over to x coordinates with the help of 


(3) 
dx, (1 x? ) 
dx 16k? 


which integrates to 





2 


x* \} x 
x= ix(1 -—) +2k arcsin — (23) 
16k? 4k 


TABLE I. Approximations of the 2s eigenfunction. 











2ar x F (29) F( 2) F 
0 0.000 0.000 0.000 0.000 
1 2.000 1.154 1.168 1.168 
4 4.000 — 0.264 0.010 0.000 
6 4.899 — 1.540 — 1.308 — 1,337 
8 5.657 — 1.861 — 2.181 — 2.164 
9 6.000 — 1.660 — 2.414 — 2.370 
12 6.928 — 0.180 — 2.642 — 2.392 
16 “8.000 1.876 — 2 — 1.758 
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The behavior of x,/k as a function of x/k is 
shown in Fig. 5. At first, they are very similar, 
but as the latter approaches 4, the former falls 
back and approaches 7 with zero slope. 

This correction brings the nodes quite well to 
coincidence. One finds now the main discrepancy 
in the size of the loops, the true loops increasing 
somewhat more sharply with r than the approxi- 
mate ones. To correct this, we introduce as our 
second step an unknown factor f into the wave 
function, writing 


F=f-, 
= x,Z1(x;). 


It is clear from the purpose of the correction 
that f may be assumed to vary slowly compared 
to ®. Substituting into (2) and using the equation 
for ®, viz., 


where 


——+-- 
Rk? 6 


dr® dr,/dr dr 
we find 


ad’ d*r,/dr*® db a® 2a 
( Jo=o, 


d*f df d’*r,/dr*\d® 
“+ | ———-}--- =0) 
dr? dr dr,/dr jdr 


According to our assumption, the first term can 
be neglected compared to the others and we get 








f ( =) (24) 
" dr gale 
where we have introduced the abbreviation 
dx; x? } 
nt=—ra=a(1- ) ; (25) 
dx 16k? 


Thus we get finally for F 
F=xx,'x,*-| J;(x,) cos kr +N,(x,) sinkw}. (26) 


It is perhaps worth while to explain the correc- 
tions of this section by physical reasoning. 
The function 

xZ1(x) 


is the correct solution of (2) for zero energy. 
According to the Wentzel-Kramers-Brillouin 
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0 . 2 3 


Fic. 5. Plotting the scale-corrected argument x, against 
the true x [Eq. (23)]. x, must be used instead of x as 
the argument in the Bessel or Neumann function to make 
our approximation quantitatively correct. It is an essential 
feature of the method that it does not extend beyond 
x=A4k. 


approximation, the true wave function for a 
negative energy state must differ from it on two 
counts: 

(a) The wave-length, obeying De Broglie’s 
relation, must be longer. This reasoning leads to 
a change in the abscissa scale which agrees 
exactly with (23). 

(b) The amplitude of the wave function is 
inversely proportional to the square root of the 
classical velocity. Hence we should multiply our 
wave function with a factor f which equals 





(— velocity at false aod) 
right velocity at right position 


This formula reproduces exactly (24). 

It is an inevitable consequence of such 
reasoning that the corrected wave function 
breaks down at x=4k because this is the limit 
of the classical orbit. But until the neighborhood 
of that point is reached, the approximation is 
greatly improved numerically. To illustrate this 
fact, the Schrédinger wave function for the 2s 
state is given in Table I, and the approximations 
(20) and (26) are compared with it. 

Fitting the wave function (26) to the wave 
function (15) at the point x = x» is now a straight- 
forward matter. One finds 


Jo(xx) sin $x4* — Ji (xx) cos $x,* +4 { (x?/x,**) — (1/24) } Ji(xx) sin 4x,* 





tan kr = — 


No(xx) sin $x,.*— Ni(x,) cos $x4* + 4 { (x?/2y**) — (1/2x,)} Ni(x,) sin hx,* 


(27) 
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x, and x;,* are derived from the cut-off coordinate 
x=x»9 by the formulas (23) and (25). xo itself is 
derived from ro through (15). 

Formula (27) is of course less convenient to 
use than (22) because the unknown is contained 
implicitly in the right-hand expression. That 
expression varies, however, much more slowly 
with & than the left (except in the excluded case 
when x is very close to 4k) and hence the formula 
is suitable for successive approximations, starting 
with an approximate k on the right. One or 
two steps are amply sufficient to get the answer. 

Applying the formula to the case of Tibbs? 
we get 
k=1.153 for the 1s state 
k=2.159 for the 2s state 


in complete agreement with his results. 


4. SOLUTION FOR OTHER STATES 


The method of Section 3 is applicable in the 
same form to states of non-vanishing angular 
momentum. We shall restrict ourselves therefore 
to the discussion of new features and the quota- 
tion of results. 

Following the same reasoning as previously 
we are led to an equation for the scale-corrected 








re—a\} (r—b)(c—r)\! 
(r—a)!—a! arctan ( ) (——) + 
a b+c 


























Fic. 6. Radial potential for a particle having non- 
vanishing angular momentum which moves under the 
influence of a Coulomb force. 
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2a I(I+1) a? 
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We notice that this expression is considerably 
more involved than in the previous case. This 
is explained with the help of Fig. 6, which shows 
the potential for radial motion. In the actual 
negative energy state the classical orbit has now 
an inner limit 5 as well as an outer limit ¢ which 
are the two roots of the denominator above; 
and the zero-energy orbit has an inner limit a, 
which is the root of the numerator. It follows 
from this that r,; will be a complex function of 
r unless the constant of integration is picked in 
such a way that when 


r=b 


then 
Tr =a. 


Such a choice has the additional advantage that 
it brings the first inflexion points of the approxi- 
mate and the true F automatically to coinci- 
dence. We find thus 









(1 1)! 
r—b\} b r 
3(b+c)! arctan (—) —a'arctan |——| . (28) 
c—r 1 1 
ir) 








The new formula is seen to go over into (23) in 
the special case when a=)=0. 

No special difficulty arises in the determination 
of the factor f. If we define as in (3) 


(29) 





x= 2(2ar,:)! 


and as in (25) 











att) xt} 
d. 2 16k? 
eed : » (30) 
dx , 41(/+1) 
q x,? 





then we get for the wave function F: 


§ «-i 
F=xxyixe: Zoiyi(Xei), 
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where Z2:,:(x) is a combination of Bessel and 
Neumann functions of order 2/+1. It seems 
logical to follow the procedure that led to (26), 
that is, to substitute for Z the linear combination 
(20). However, such a procedure is incorrect be- 
cause the substitution of x,; for x has brought in 
some constant factors. x,:, through (28), (29), 
and (3), is a regular function of x right through 
the inner limit point of the orbit x,.=2[/(/+1) }}; 


the same applies to xf, as defined by (30). They 
both are also equal to zero simultaneously with x, 


: > i ‘ ee 
and their ratio x;:/x,2, tends to 1 in that limit; 
but the ratio x;./x does not do so. We find 


[1(+1) }! { A/ (00-41) 19) 4-19 
awe 
© 
[1(L-+1) Jy t/t 
(:-———) - 
; k 





° 








J 


which is a quantity lying between 1 and 0.8578. 
Now the Bessel-like and Neumann-like functions 
are characterized by their behavior at the origin 
where they are to behave as 


(3x)?*1 (21) (3x) 
—— and —- 
(2/+1)! 7 





’ 





(k+))! 
(k—1—1)!k2+1 





tan kr= —C;,, 
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TABLE II. Comparison of the 4d wave function with 
approximation (32a) 
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respectively. Hence, if we define 
| (1(d+-1) Pry ery 
1-——) | 
( k 


[l(l+ 1)} { e/(0C1+1))4y 41 
‘omic ie, 


9 


e- 





Ciu= 








we get F equal to 


_ + «-—} ; 
F = XXEIXK1 Cord or4-1(Xe1) cos kr 


_y(R—1—1) te24+! : 
Cx (e+)! Noi41(Xe1) sin bet. (32a) 
A comparison of the Schrédinger wave function 
with its approximation (32a) is carried out in 
Table II. The factor C,; equals- 1.424 in this 
case and could certainly not be neglected. 
To get our equation for k we must fit (32) to 
(14) at the cut-off. The result is 





x.* 
J 21( 1) J i44($xu*) — Fars sn1)Jra($xu*)—+ WJ anya ur) Tess (hxe®) 
x 


kl 





(33a) 


* 


; XK 
Nai(xnr) Ji +4(3x4*) es Noug.1(2¢1) T1442") —+WNo141(et) J sxe") 
x 


The factor ¥ equals 


2i+1 2/+1 1 


+ igs 
x x 2(1-—= 
x? 





kl 


1 x? x,*? 
x (—+5,-* 


x ££ x*3 
kl kl kl 


(34) 


In spite of appearances to the contrary, the 


kl 





factor ¥ tends to zero for small x and remains 
finite for x,=2[1(/+1)]!. The variables entering 
into (33a) and (34) are defined in (3), (25), (28), 
(29), (30), and (31). 

It is an interesting accident that the factor 
entering (33a), namely, 


(k+I)! 
“"(k—1—1) Wet)" 


is almost 1 for all possible values of & and /. 
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For instance, for k=2, /=1, we find 1.011, for for F, apart from a constant factor 
k=3, l=1, we find 1.002 and so on. It reaches is 

its largest value for k=, J=k—1, namely, *=**&ee {J2r41(%e1) Cos ka 
1.074, while for =k —2, it is smaller than 1.026. + Noei4i(xx:) sin ke} (32b) 
Thus we may often neglect it entirely, and write and instead of (33a) 











x * 
So1(%41)Fr44(30%*) — Tata (nr) es bee) + WSa043 1) Tea hoe") 
kl 


tan kr= — ~ : (33b) 
x 
Noi(xx1) J 144(3%%*) — Navex ui) Jia Qee*) + Wane a(a) +a(3x2*) 


kl 

















If we apply Formula (33) to the 2p state It seems to the author that the elimination of 






studied by Tibbs? we find purely numerical work in quantum mechanics is 
very desirable because of the greater generality 
k= 2.004. gained. The present paper makes the cut-off 







Coulomb field an analytically soluble case except 

for energy levels which are so close to the cut-off 

k=2.005. in the potential energy that only an exponential 

tail is situated beyond the cut-off radius. More 

Unfortunately a check with his result does not precise criteria, together with formulas to replace 

mean very much here because his value may lie (21), (27), and (33) will be given in another 
anywhere between 2.00 and 2.041. paper. 





The approximate formula (21) gives 
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According to the phenomenological theory of quantum fields, fluctuation phenomena corre- 
spond even to the lowest energy state of a system. The average value of field quantities is zero 
in the case of a vacuum. This value changes, however, if we introduce a particle into the field. 
The determination of the average value of the field fluctuation leads immediately to the expres- 
sions of the proper field of an electron. The calculation becomes considerably simplified by an 
appropriate use of Dirac’s theory. 













A. INTRODUCTION an electron. Even the value of kinetic energy of 
an electron cannot be obtained from classical 

1. Purpose of Paper vai : ; 
50 of Gin Page physics in a satisfactory way since the value 





T is well known that classical theory does not given by the mechanical picture does not agree 

represent, so far, a consistent scheme of at all with the value one would expect from the 

physical phenomena, even if we apply it to a_ transformation properties of static electromag- 
simple problem as the problem of the motion of _ netic fields. 
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THEORY OF STATIC FIELDS 


The first part of this paper intends to study 
the proper field of the electron from the point of 
view of quantum theory, in order to see whether 
or not quantum theory leads to modifications of 
the transformation properties of microscopic 
static fields which are likely to attenuate the 
difficulties met in classical physics. It answers 
at the same time the question in which way 
quantum theory accounts for static fields which 
are responsible for elementary macroscopic 
phenomena as electro- and magnetostatics. 

The paper will give an account of a number of 
investigations which have been carried out 
during the last years by the author and a few of 
his collaborators and which could not yet be 
published in-detail because of the events in 
Europe. 


2. A Simple Example 


In order to have a general view of the prop- 
erties of the formalism of quantum field theory, 
we consider first a scalar field, described by a 
complete set of orthogonal functions in an 
arbitrary large cube of volume L’*, 


uy, =e**/L). (2.1) 


Px = ey ; 
These obey the Klein-Gordon equation 


1 0°¢ 
Ag-—— ——_ «*9=0. 


(2.2) 
c*? df 


The static solution of Eq. (2.2) is given by 


g= —Ae*'/r. (2.3) 
We shall see how quantum field theory accounts 
for this solution. 

According to the formalism of quantum field 
theory, we have to consider every single Fourier 
component (2.1) as an oscillator, the energy of 
which can take only integral multiple values of 
he,/2x. If no particle is present in the field, 
the Hamiltonian of the field is given by 


H= Sx (pu? +0r2gu?—hox/2e), (2.4) 


where p, and gy are canonical variables, 


PuQu — Ve Pu =h/ 221. (2.5) 
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Then, the field is represented by the operator 


e=2 i(tn)'¢| (a+r) 


1 
+(a-—m)u|. (2.6) 
WW 

It is well known that the Hamiltonian (2.4) 
describes fluctuating variables p, and qx, the 
average values of which are equal to zero. The 
average value of the field (2.6) in free space is, 
therefore, equal to zero. 

Introducing an infinitely heavy particle with 
the coordinates R into the field, we can consider 
the kinetic energy of the particle to be a constant 
and the Hamiltonian of the system can be written 


H= | 3 (Pu? + wn?gu? — hax / 22) 


1 
+3(4)'c-A | (a+—») -ux(R) 


IW 


(ota) we] 


= Die 3 (Put Pu”)? +on*(Get+ Qn’)? 
— heo,/ 2 — px? — wr2gn}. 


(2.7) 


The constant A plays the role of an interaction 
constant or “‘charge’’ in this theory. The ab- 
breviations p,° and q,° are determined by 
ux (R) —u,*(R) 
p,°= (43r)'c-A ° . : 
21, 





ux(R) + u,*(R) 


2 





qx° = (4) ic-A ° 


2a% 


The Hamiltonian (2.7) represents an oscillator 
system, the average values of its canonical 
variables being shifted by the (negative) amounts 
given by (2.8). Introducing these average values 
in the expression (2.6), we obtain immediately 
the average value of the field 


ux*(R)ux ux(R)ux* 
+ 


2 2 





@=—2nc*?A 


k Wk Wk 


e~*it—Rl 
=-—A (2.9) 


lr—R| 
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in agreement with expression (2.3), since 
(w/c)? =k? +x? 
de esik(r-R)  g-«lr—RI 


>» = 
rR+e2 L? |r—RI 


and 





The two last terms in (2.7) represent the 
infinite self-energy of our field, which, however, 
is without any interest for our purpose. 

Replacing our simplified field model by the 
oscillator system representing the electromag- 
netic radiation field and the particle by a Dirac 
electron, the present scheme permits immediate 
calculation of the proper field of the electron to 
the approximation in which the electron is 
sufficiently heavy not to be appreciably in- 
fluenced by the field. 

Before, however, entering into the discussion 
of the proper field of the electron, we shall 
develop a mathematical apparatus, which will 
not only be helpful in the following calculations, 
but also simplify the considerations which we 
will have to present in the later parts of this 
series of investigations. 


B. MATHEMATICAL APPARATUS 
3. The Solution of Dirac’s Equation 


We shall write Dirac’s equation in the fol- 
lowing form: 


1ayv 1B 
——-—-+ea grad ¥+—V=0. (3.1) 
c Ot A 


The constant A can be expressed by 
A=h/2xmc, (3.2) 


the matrices a and 8 are supposed to have the 
usual form. 

We shall call a solution of Eq. (3.1) a complete 
set of orthogonal functions, depending on the 
variables x, y, 2, t, a, and 8. One verifies easily, 


that 
pot+me— yp 
{ L32po(po-+me) } } 


2ni 
“exp |= rnet-+r) (3.3) 





is a complete solution of Eq. (3.1) if p denotes 
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an arbitrary vector in space and 
po? — p? =m?*c*, (3.4) 
7 =Ba. (3.5) 


We shall refer to the form (3.3) of the solution 
as the standard form. 
In Table I the solution (3.3) is written ex- 
plicitly, without the common normalizing factor. 
Every column of (3.3) represents a set of four 
spinors which is generally called a partial solution 
of Dirac’s equation. The four partial solutions 
included in (3.3) correspond to two sets of 
plain waves moving in opposite directions. This 
is an essential feature of the standard solution. 
While the asterisk (*) denotes in the usual” 
way the conjugate complex of a quantity, we 
shall denote by the sign “‘til’’ (~) the conjugate 
transposed of a matrix. The conjugate trans- 
posed of Eq. (3.1) is 
law oe 
——-—+grad Va——V6=0, (3.6) 
c at A 


2m 


-—"apuct+or) 
h 


Vp=exp 


po+mc+yp 


x _ (3.7) 
{ L32po(po+mce) }! 





One verifies easily the normalization conditions 


[Fe ¥dr= sy, (3.8) 
a 1 

hal ar (3.9) 
" 1 


Dirac’s matrices describe quadrivalent vari- 
ables, the four values of which we shall consider 
to belong to different values of the character of 
the electron. The character includes two bivalent 
sets of values, (a5) and (6y), belonging to 
different spins, and two bivalent sets, (a8) and 
(v6), belonging to different signs of the electron. 

An electromagnetic field splits, for an observer 
in a given system of reference, into two parts: 
the electrit and the magnetic field. In a some- 
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TABLE I. Solution of Eq. (3.3). 











2xi 2ni 
exp [=e ~pu | exp [ x — pot) 


2ni 2ni 
| exp [Fee + out | exp Fe +o 





a Pot mc 0 
B 0 Potme 
Y —Pps = (pe—tpy) 
6 = (pet+ipy) +p: 
a B 





+p: + (P:—ipy) 
+ (p2+1p,) —?p: 
Pot mc 0 
0 Pot mec 
Y 6 








what analogous way, the character of an electron 
splits into spin and sign. 

The relation (3.8) accounts for the orthogo- 
nality in space, the relation (3.9) expresses the 
orthogonality in character. 


4. Dirac’s Matrices 


It is well known that Dirac’s matrices 1, 6, @ 
form a fundamental system of matrices. They 
permit the establishment of a complete linear 
basis system of 16 matrices which permits the 
expression of any matrix of four rows and 
columns by a linear expression. We shall denote 
the linear basis system as follows 


B 
1 —«a 
17 — Be (4.1) 
T —o¢ 
iBT 


with the following notations 
i 
y=6a; o= ~sreategi t=4(a-o@). (4.2) 


The matrices (4.1) are of the Hermitian type. 
We shall call matrices of the type 


Oo O Oo O 
Oo O Oo O 
Oo O O O 
Oo O Oo O 


respectively matrices of first and second kind. 
The matrices 


1, 8B, @, Bo 
are of first kind; 
a, ty, 7, tBr 


of second kind. 





In general, a matrix will contain parts of both 


kinds. We shall denote by 


\M=3(M+ MB), 
2M =}(M—BMB) 


(4.3) 


the part of first and of second kind of the 
matrix M.! 


5. The Cartan Transformation’ 


We introduce in Dirac’s equation a general 
interaction term of Hermitian type, U= U 


1 dV 1 2m1 
—-—+a grad ¥+-6¥+—UwW=0. (5.1) 
c Oot A h 


The conjugate complex equation to (5.1) is 


1 oy* 1 21 
———+a* grad y* ——6¥* ——U* ¥* =0. (5.2) 
c at A h 


Introducing Cartan’s matrix 


+1 
p= —ifBa,= - =" , pat, (5.3) 


+1 


multiplying (5.2) from both sides by p and 
putting 
b= p¥*p, (5.4) 


we obtain, because 
pBp=—B, pa*p=a, 


1A more detailed investigation on the properties of 
Dirac’s matrices will be found in Fernandes de Sa’s thesis, 
Pérto 1943, in which the general properties of the micro- 
scopic space-time continuum are discussed. We shall 
have to refer frequently to Fernandes de Sa’s paper. 

? The author is indebted to Mr. Jean Pirenne, Lyon, 
for having drawn his attention to this important transfor- 
mation, given by E. Cartan. 
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the equation 









1 0d 4 
—-—-+e grad $+-— 
c Ot A 


2a 
iat anne (5.5) 
In the general case, we have to put 
U=V+W, 
V=Ao+eA+iyE — eH, (5.6) 
W=6B+7rC,—eC+i8rD, 


where the A, B, C, D, E, and Hi’s are real 
quantities, forming appropriately chosen tensors 


of a Lorentz transformation. Since 
pV*p=V; pW*p=—W 


Eq. (5.5) takes the final form 





1 0@ 1 
—-—-+a grad $+-—66— 
c at A 


2x1 
—(V—W)#=0. (5.7) 
h 
Cartan’s transformation provides us with a 
convenient means to study, with the help of 
Eq. (5.7) the behavior of states, which, according 
to the usually considered Eq. (5.1), belong to 
negative energy values. In the restricted case, 
in which only an electromagnetic field, repre- 
sented by the four vector Ao, A is considered, 
we have to postulate that no additional assump- 
tion be introduced, which destroys the necessary 
symmetry corresponding to Cartan’s transfor- 
mation (Dirac’s “filling up rule’’). 








6. Spatial Rotation and Lorentz 
Transformation 


Since we will not have to use these transfor- 
mations explicitly, only a brief survey shall be 
given here. The general theory is contained in 
Fernandes de Sa’s quoted paper. 

A spatial rotation can be represented in this 
formalism by a left-hand operation 






W’=RV. (6.1) 





The transformation matrix R is of first kind; 
it depends on Eulers angles ¢, 3, ¥ and on the 
matrices 1 and ¢, R(¢, 3, ¥; 1, @). 

An essential new feature appears as soon as 
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we consider Lorentz transformations, 












W’=L,Y, (6.2) 
where 
“uv u 
L,= Ch-———aSh-, (6.3) 
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represents an ordinary spinor transformation. 

We have to note that the spinor transforma- 
tion (6.2) destroys the standard form of our 
solution. We have, therefore, after a spinor 
transformation, to rearrange our solution by a 
unitary transformation in order to reestablish 
the standard form. Physically, this fact means 
that the transformation properties of a fluctu- 
ating quantity of second kind does not obey the 
ordinary transformation laws under a Lorentz 
transformation. 



























7. General Motion 






We shall call a general motion a complete 
set of orthogonal functions 





Sees mite f bpbdr. (7.1) 





The transformation coefficients are matrices of 
four rows and columns, and depend explicitly 
on time and form a unitary transformation. If 
we decide to carry out summation on double 
indices without mentioning explicitly the sign 
of summation, they satisfy the conditions 








Sip'SPe= Sie} SPs Sip = Sprp- (7.2) 
While in (7.1) Vp represents a solution of 
Dirac’s Eq. (3.1), ®; does not satisfy this equa- 
tion. In general, it will be possible, however, to 
find a function U depending explicitly on time 
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in order to have ® satisfying an equation of the 
type of (5.1). 


C. THE PROPER FIELD OF AN ELECTRON 


8. The Radiation Field 


In order to account for the existence of 
photons, which according to the phenomeno- 
logical theory has not to be explained but has 
simply to be taken as a postulate, we introduce 
the operators 


1 @ 
Ox, Ny=-— (8.1) 
4 00% 


operating on eigenfunctions of the type 
, 1 
r =exp i £ (ma v; In 2), (8.2) 
” 


where the m,’ can have any positive or negative 
integer values including zero. If we want to 
interpret the numbers m,’ as numbers of photons, 
we have to exclude all solutions (8.2) in which 
all the values of m,’ do not have the same sign. 
We shall denote by m, positive integer values 
including zero. 
If we denote by 


We =hw,/ 2a (8.3) 


the energy of a photon corresponding to a wave 
vector k, the eigenvalue problem of the radiation 
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field becomes 


Oo 2m 
(+-+— Dix wee )r =0. (8.4) 
at oh 


We have to attribute to the field energy the 
same sign as to the energy of the particles we 
want the field to interact with. We have, there- 
fore, to introduce a scalar operator ¢ which 
represents the sign of the electron charge, &=1. 
We shall suppose « to commute with Dirac’s 
matrices. It is, for our considerations, an essential 
fact that Dirac’s theory does not contain any 
operator of the type required here. Phenomeno- 
logical theory has, therefore, to introduce this 
operator ad hoc, in an unsatisfactory way. 

The solution of (8.4) may now be written 


2x 1 
Il'=exp : | meret—emde + In an} | (8.5) 
k 


The vector potential of the field will be 
determined below. We note, however, that the 
radiation field is described in classical theory by 
a complete set of orthogonal functions 


he 





Ayu, (8.6) 


A, = 
THE 
where A, is a unit vector perpendicular to k, 
determining the polarization of the plane wave. 
ux has already been defined by (2.1). The 
amplitude in (8.6) is normalized to correspond 
to the energy of one photon w,. 





9. The Interaction Energy 


For the interaction energy we can use the usual expression, but we need to take care that only 
electrons and photons of the same energy sign interact with each other. The interaction energy 


becomes 


he 
U(e) => —e- Aya: 
k (2 





WWE i 


{3(1—«)-[(Vx)!-exp (70%) -ux*+exp (—70x) -(Nx)*- ux] 


+4(1+¢e) ’ [(N,*)!-exp (—16x) ‘Uy +exp (76x) . (Ny*))-u,*]}. (9.1) 
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which is defined by 


The operators \/N and \/N* can be expressed by a differentiation of the order } 
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With the help of (9.2) the interaction energy (9.1) can be written 








7 
exp (7 
he 
U(e)=>0—e- Aya: . {Dy} exp [ie(kr— 6.) ]+exp [—ie(kr—6,) ]Dx4}. (9.3) 
. (2aw,)! L} 
and we can consider 
T 
exp (i) 
he 4 , 
A=> A {D,} exp [ie(kr— 6.) ]+exp [ —ie(kr—6,) |D,'} (9.4) 


. Gan 2 


to be the vector potential operator of the field. 


The wave equation of the complete system becomes 


| re) 18 2x 2n1 = 
——+ca grad+—+— U(e) -—— > wn: My |¥=0 (9.5) 
| ot A h hk 
and transforms under a Cartan transformation to 
0 18 2m 21 _ 
——+ce grad + ——U(—) ——™ F wn Ns] =0. (9.6) 
| ot A h h * 


The functions VW and @ depend not only on the 
electron variables x, y, 2, t, a, 8, but also on the 
field variables ¢, 0,. The expressions (9.5) and 
(9.6) account correctly for the relation between 
electron and positron, while the usual interaction 
term, together with Dirac’s filling up rule, does 
not account in a correct way for the proper field 
of an electron. 

It can be objected, however, that our inter- 
action term (9.3) is not in agreement with 
experiments since it does not contain any rule 
which permits the exclusion of spontaneous 
radiation transitions between states of different 
energy sign. This objection cannot be main- 
tained. The quoted transitions occur only in 
second order while we shall show here that (9.3) 
applies only to non-relativistic first-order prob- 
lems and that not only the interaction term, but 
even the whole scheme of quantum field theory 
has to be abandoned. The question of the quoted 
radiative transitions does, therefore, not appear 
explicitly in our considerations. 


10. The Field Perturbation 


We have now to solve, by an appropriate 
method of approximation, Eq. (9.5), and we 
must determine the average value of the ground 





state of the field (9. 4). Supposing that no free 
photons are present in the field and that the 
electron has the momentum p, we can represent 
the solution of (9.5) in first approximation by a 
series which corresponds to states in which one 
“virtual” photon has been emitted, 


V=W,-exp (—E } In 27)+Dd Vp 
p’k 


21 
“exp (~ east —ieth— 4 In 2m Jer (10.1) 


where the coefficient matrices c,, remain to be 


determined. 
Introducing (10.1) into (9.5) we obtain 


immediately 
dc wth 'k 


2ri ; . 
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ew exp (— > 3 In 27) 
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and, with the help of relation (3.8) 
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In order to evaluate the relation (10.2), we our restriction limits our calculation of the proper 
remember that we are interested only in small field of the electron to large distances of the 












































values of momentum particle, r>A. For the present calculation, it 
h implies, however, a considerable simplification, 
) P, P, 5 e<me. (10.3) since we can put 
T 
, — bo= po’ =mc 
rhe restriction (10.3) is due to the fact that we 
cannot expect our method to give reliable and since we can neglect all terms which are of 
results as soon as the particle cannot be con- higher than first order in the momenta indicated 
) sidered to be very heavy. We shall see later that in (10.3). One verifies easily that 
dcp 2a 2m 
oS <a * fx p (-—“eaut) by aW pe'***dr 
dt h (2aw,)! L} 
2xi  ehc Ax 2mi yp’ YP 2x1 
) = —— exp | ~—"(ome:+-aa)s}| (14+ )a(1-)} exp (amc) 
h (2mw,)! L! h 2mc 2mc h 
2ni ehe *| 2ni B 7 . 21 
= ——— ——— —}{ € exp —29met) ———|[p’+p+i(eXp’—p) |] exp (-—"«w) . (10.4) 
h (2mw,)' L} ( h 2mc h 
with 
) h 
p—p’+e—k=0 
‘ 2r 
2 and, finally, 
ehc {Aya-B [2mi 3 h 2ni 
Cox = — | exp [ase —ean)| A,| 29+i( oxe—t) | exp ( ~~ wt) . 
1 (2aw,L*)'| 2mc* h 2mcew), 2x h 
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11. The Average Value of the Field 






We now have to calculate the average value of the field (9.4), 







- — 
—A= | VAVd9. (11.1) 
L3 









We observe that the field A depends on coordinates r in space while the functions ¥ depend on the 

electron coordinates, which in this paragraph shall be denoted by capital letters R. The average 

field depends, therefore, on time, on the coordinates r and R, and on the variables ¢, a, and 8. 
Introducing in (11.1) the respective values from (9.4) and (10.1), we find 
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In our approximation, the functions V can be replaced, a part from a factor L~!, directly by the 
respective exponentials, 
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where we introduce the values of the coefficients c,,, found in (10.5). The result of the calculation is 


A=A.+A,+A. (11.2) 
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12. The Summation of the Series 


The series (11.3), (11.4), and (11.5) can be evaluated. We note that the summation has to be 
carried out not only on the wave vector k, but also on both directions of polarization. One veri- 


fies easily the following relations 
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With the help of (12.1) and (12.2), we find immediately 
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13. The Proper Field of the Electron 


We now have to interpret the expressions 
(12.3), (12.4), and (12.5). Since the velocity of 
an electron corresponding to our standard solu- 
tion (3.3) is equal to —v, A. corresponds 
directly to the magnetic field because of the 
translation of the Coulomb field 

e 
Metequncmme, (13.1) 
|r—R| 
The gradient in the vector potential (12.3) does 
not contribute to the magnetic field and has no 
physical signification. 
(oXr—R) 


PEN avira ad 13.2 
eup TOR) (13.2) 


a(r—R) 2mi 
— grad —— | exp (—29met), (12.5) 
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b= 
4amc 





represents the vector potential of the magnetic 
spin field of the electron. 

The last field vanishes, since 8 anticommutes 
with @ 


A. =0. (13.3) 


The proper field of the electron consists, there- 
fore, in our approximation of two parts which 
belong to the charge and to the magnetic moment 
of the electron. It obeys Maxwell's equations 
and accounts correctly for electrostatic and 
magnetostatic phenomena. 

The field A, is of second kind. It corresponds 
to excitations of the field due to virtual pair 
production. Its average value vanishes. Still, 
we know that this field plays an important part. 
Mr. Jean Pirenne has made a thorough investi- 
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gation about the way in which the field A, 
contributes to the interaction of electrons. He 
finds, in agreement with (13.3), that no contri- 


butions to interaction result at finite distances. . 


Still, the field A, gives rise to interaction terms 
between electrons and positrons at very small 
distances. These interaction terms are of the 
order of spin-spin interaction and can be inter- 
preted as forces of the Heisenberg and Majorana 
type between unlike particles.* 


14. Thibaud’s Experiment 


The next task is to investigate the terms of 
higher order of the electron field. The higher 
order terms which have been calculated‘ are 
rather complicated and no longer obey Maxwell's 
equations. If we want to attribute a physical 
meaning to them, electrodynamics becomes an 
open scheme which does not longer permit any 
theoretical predictions. If the difficulties con- 
cerning the transformation properties of static 
fields are due to deviations of Maxwell's theory, 
the question arises whether or not Maxwell's 
theory represents a complete description of the 
electromagnetic field in vacuum. 

Owing to Professor Jean Thibaud’s kind 
interest in the problem, it has been possible to 
examine this question from the experimental 
point of view. Let us suppose that a. small 
fraction of the electromagnetic field is not 
described by a vector potential but by an 
antisymmetric tensor potential, as some of our 
quoted calculations indicated. It can be easily 
shown that such a potential corresponds to the 
existence of longitudinal photons of very high 
penetrating power which could, however, be 
detected if their wave-length approaches Comp- 
ton’s wave-length. Professor Thibaud examined 
the radiation due to electron impacts under a 
voltage of about 280 kv. The penetrating power 
of the radiation has been studied with a thick 
layer of lead, up to 140 mm. The result was nega- 
tive; no indications for longitudinal photons were 
found; the penetrating part of the radiation corre- 


+ Jean Pirenne, Thése, Lyon (1941). The author has no 
og om of knowing at present whether or not this paper 
as been published in the meantime. See, however, 
Cahiers de physique 4, 1, 6 (1941). 
*One part of these calculations has been made by 
Mr. S. Héxter. 
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sponded exactly to the known cosmic radiation.® 

We concluded from Thibaud’s results that the 
origin of the fundamental difficulties concerning 
the transformation properties of static fields is 
not the electromagnetic part of the problem and 
has to be searched for elsewhere. 


15. Conclusions 


The following conclusions can be drawn from 
the first part of this paper: 

1. Elementary electrostatic facts lead to a 
consistent mathematical scheme which can be 
applied to Dirac’s theory. 

2. The phenomenological scheme of quantum 
field theory is not appropriate to our problem 
and has to be abandoned. 

3. The proper field of the electron depends 
essentially on the character of the electron. 

4. No progress of the theory is possible with- 
out a careful study of the kinematical behavior of 
the electron character according to Dirac’s theory. 

That study shall be communicated in the 
second part of this paper. 


ACKNOWLEDGMENTS 


The first part of this paper was completed in 
France. The author desires to express his re- 
spectful thanks to Professor Paul Langevin 
(Paris) for the possibility of working in Paris in 
1938. His thanks are due to his friend Jacques 
Salomon (Paris). The author thanks respectfully 
Professor Louis de Broglie for the support given 
him. The author is very much indebted to 
Professor Jean Thibaud (Lyon) for. his kind 
interest in this work, and it is a great pleasure 
to him to express his gratitude to his friend, 
Professor Max Morand (Lyon), who did not 
miss one single opportunity, official or private, 
to facilitate this work. The author desires to 
thank his young friends, Jean Pirenne (Lyon) 
and Fernandes de Sa (Porto) for the enthusiasm 
with which they assisted this work under, 
indeed, very difficult conditions. Last but not 
least, the author expresses his gratitude to Pro- 
fessor Enrique Gaviola (Cérdoba) for the possi- 
bility of continuing this work at his observatory. 


5 The author is indebted to Professor Jean Thibaud for 
the communication of his unpublished results by letter of 
November 7, 1940 and for his kind permission to refer to 
his measurements. 












PHYSICAL REVIEW 


Letters to the Editor 








ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 





New Planetary Systems 


H. K. Sen 
Allahabad University, Allahabad, India 
October 1, 1943 


BSERVATION has at last succeeded in giving us a 

verdict on the vexed problem which had been engag- 
ing the attention of the theorists for some time—whether 
our solar system is the only one of its kind in this vast 
universe. From the study of parallax observations Dr. 
K. Aa. Strand! has recently been able to prove the exist- 
ence of a planetary body companion to the 61 Cygni sys- 
tem. The system shows a periodic deviation from the 
Keplerian orbit, and this can only be explained on the 
supposition that a third body (61 Cygni C) of small 
mass, about sixteen times that of Jupiter, revolves around 
one of the stars of the binary system. The luminosity 
of 61 Cygni C is found to be very low, and it can very 
well be classified as a planet. 

A parallel discovery to 61 Cygni C has also been almost 
concurrently made. Reuyl and Holmberg* have discovered 
a similar planetary body of even smaller mass in the binary 
star system, 70 Ophiuchi. It is worthy of note that neither 
of the stars 61 Cygni and 70 Ophiuchi exceeds a distance of 
10 parsecs from the sun. 

We have, at last, evidence of the discovery of two 
planetary systems other than our own, and, what is more, 
these ‘happen to be very near our own system. The sugges- 
tion naturally comes to mind that we may by similar 
means be able to find many more planetary systems in our 
galactic system, and the number of planetary systems in 
this universe may be quite large. It will be interesting in 
this connection to review briefly what theory has to say 
on this question. 

On the tidal theories of Jeans* and Jeffreys,‘ planetary 
systems should indeed be rare. Taking the generally 
accepted age of the universe as of the order of 10° to 10” 
years, and Jeans’ estimate (1929) of the frequency of 
planetary systems on his theory—about one per 5000 
million years—we have at the most two planetary systems 
in our galactic system. On Lyttleton’s® binary star collision 
theory, the probability of the formation of planetary sys- 
tems is almost nil. 

The probability is much greater on Banerji’s* Cepheid 
theory (1942) in which a passing star is supposed to break 
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up an oscillating Cepheid. Encounter at only a moderate 
distance is necessary in this theory—it need not be grazing 
or close. Banerji has stated his conclusion in the following 
words, “One conclusion seems to be irresistible. If the 
theory be correct in its essentials, there may be more 
planetary systems than at present supposed.’’ The obser- 
vations thus definitely favor Banerji’s conclusion. 

It may be remarked that the presence of planetary 
systems in binary stars is very plausible on the Cepheid 
theory of the origin of binary stars recently advanced by 
the author.’ A Cepheid is supposed to break up into a 
binary star system on account of increase of angular 
velocity due to radiation of energy. The filament connect- 
ing the two stars affords a chance for planets to be formed 
by condensation. 

'K. Aa. Strand, Pub. Astronom. Soc. Pac. 55, 29 (1943). 

2? Reuyl and Holmberg, Astrophys. J. 97, 41 (1943). 

3 J. H. Jeans, The Problems of Cosmogony and Stellar Dynamics (1919), 
p. 275; Astronomy and Cosmogony (1929), pp. 7, 22, 409. 

4H. Jeffreys, The Earth (1929), p. 16. 

5 R. A. Lyttleton, M.N.R.A.S. 98, 536 (1938). 


* A. C. Banerji, Proc. Nat. Inst. Sci. Ind. 8, 173-197 (1942). 
7H. K. Sen, Sci. and Cult. 7, 582 (1941-1942). 





Quantized Probability 


ALFRED N. GOLDSMITH 
580 Fifth Avenue, New York, New York 
November 13, 1943 


HIS communication proposes for development the 

theory that there exists a finite unit of probability, 
representing the lower limit of possibility. (A suggested 
designation for this unit is the ‘“chance-quantum,”’ abbre- 
viated cq.) Thus, if the probability of an event is equal to 
or greater than one cq., it may ultimately occur. If its 
probability is less than one cq., it will never occur. The 
latter statement constitutes a break with classic probability 
theory. 

Probability-quantization theory results from considera- 
tion of hypothetical experiments wherein, for example, a 
sufficiently great multiplicity of monkeys operating at 
random an adequately large number of typewriters will 
ultimately reproduce master-works of literature. Inherent 
aversion to this conclusion has led some to hold that the 
assumed result would necessarily be interpreted as fraudu- 
lently achieved. This mode of escape from the assumed 
conclusion can perhaps be avoided by postulating that the 
probability of thus reproducing great literature is so small 
as to fall below one cq. Accordingly the assumed result 
becomes flatly impossible. 

If probability-quantization is an existent condition in 
nature, Certain interesting theoretical consequences appear 
to follow. These include: 

1. If an event has a calculated probability of less than 
one cq., it will not occur. 

2. Assuming the validity of Bayes’ theorem, events can 
occur only if both their existence probability and their 
likelihood equal or exceed one cq., and if their product 
exceeds one cq. 

3. Events may differ in probability only by an integral 
number of cq. 
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4. The occurrence of events having a probability of one 
or relatively few cq. will be not only extremely infrequent 
but may have an unusual form of random distribution. 

5. For an event having an appreciable probability 
(equivalent to many cq.), a change in surrounding condi- 
tions leading to a computed change in probability of less 
than one cq. will in fact cause no change in the probability 
of the event. 

6. The interrelations between the cq. and the energy 
quantum (as well as the elsewhere proposed time quan- 
tum) may be close and significant. 

An absolute proof of the existence of the cq. faces the 
usual difficulties of proving a universal negative. At best, 
the establishment of a high probability for the existence 
of the cq. is attainable. A possible experimental method 
might be based on the following thoughts. Suppose that 
there exist » possible but very improbable events of a 
similar and identifiable type. If the selected one of these 
events not only does not occur in m trials, but also does 
not occur in mn trials (where m is a large number), there 
is a likelihood that it will not occur at all and that its 
probability lies below the assumed unit cq. 

Such an experiment might be attempted using, as a 
basis, mechanical, electrical, or electronic coincidences. 
Illustratively, there may be used a set of coaxially rotating 
disks, propelled pneumatically along their peripheries, 
mounted in low friction bearings, and mechanically isolated 
from each other. Each disk bears a circular perforation of 
the same radius, and at the same distance from the disk 
center. A collimated light beam, parallel to the rotational 
axis, can pass through the perforations when lined up, and 
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thereafter be registered by conventional photoelectric 
means. Control of perforation size, perforation distance 
from the axis, and number of disks, will enable selection 
of desired probability parameters. If a coincidence of ex- 
tremely small probability did not occur in a time which 
was a large multiple of the time of its probable occurrence, 
the hypothesis of probability-quantization would be 
strengthened. Great care might be necessary to avoid 
confusing quantum effects in the physical system with 
quantization of probability. 

If the cq. theory is found valid, interesting consequences 
may follow. Philosophically, the universe will be a system 
subject to prediction (determinism) for events of prob- 
ability less than one cq. and to non-prediction (‘‘free will” 
or normal chance) for events of probability above one cq. 
Probability theory for ordinarily likely events will be 
little affected. But the probability of small-scale and ex- 
tremely unlikely effects may be largely controlled by the 
new theory—a matter of considerable interest in certain 
modern microphysical realms. And the comparison be- 
tween continuous probability and quantized mechanics, 
on the one hand, and quantized probability and continuous 
mechanics, on the other hand, may well merit study. 

It is hoped that others will see fit to develop probability- 
quantization theory and its implications. The author par- 
ticularly wishes to express his deep appreciation of the 
searching and constructive comments on this theory which, 
in the course of correspondence, he has received from 
Professor F. Bernstein, Dr. Thornton D. Fry, and Pro- 
fessors Harold Hotelling and Richard von Mises, as well 
as other interested men of science. 
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MINUTES OF THE MEETING OF THE NEW YORK STATE SECTION 
HELD AT ITHACA, NEw YorxK, NOVEMBER 6, 1943 


HE Fall Meeting of the New York State Section of the American Physical 
Society was held at Baker Hall and Rockefeller Hall on Saturday, 
November 6, 1943. An attendance of over ninety members and visitors was 
recorded at the sessions. 
The following program of invited papers was presented: 


Morning Session 

Address of Welcome. G. H. SaBine, Vice President, Cornell University. 

Education of Physicists for Industry. J. T. LyrrLeton, Corning Glass Works. 

Metallic Mirrors by the Evaporation Process. A. FoGELSONGER, Evaporated Metal Film Cor- 
poration. 

Coulomb’s Laws, What They Are and What They Are Not. C.C. Murpock, Cornell University. 

Order and Entropy at the Absolute Zero. K. K. Darrow, Bell Telephone Laboratories. 


Afternoon Session 
Discussion: The Army and Navy Training Programs 

Army Basic Course. H. E. Howe, Cornell University. 

Navy Basic Course. P. R. GLEASON, Colgate University. 

Premeteorology Course. G. H. CAMERON, Hamilton College. 

Preinduction Training. P. G. JoHNsON, Cornell University. 
W. R. FREDRICKSON 
Secretary 
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New York, OcTOBER 23, 1943 


HE first meeting of the Metropolitan Section of the American Physical 

Society for the season 1943-1944 was held on Saturday, October 23, 1943, 

at Hunter College, New York City. About seventy-five people attended. The 
programme was as follows: 


Afternoon Session 


Contributed Papers: 

Interpretation of the Spinning Electron with Bipolar Coordinates. MARY HEWLETT PAYNE, 
Brown University. - 

Absorption of Gamma-Radiation in Terrestrial Material. C. W. SHEPPARD, Massachusetts 
Institute of Technology (By title). 
Invited Papers: 

The Study of the Ionization of the Atmosphere by Cosmic and Local Radiation. V. F. HEss, 
Fordham University. 

Production of Neutrons by Cosmic Radiation. S. A. Korrr, New York University. 

The Effect of Centrifugal Force on the Electromotive Force of Galvanic Cells. D. A. Mac- 
INNES, Rockefeller Institute for Medical Research. 

The Field Emission Arc in its Application to High Speed Radiography. C. M. SLack, Westing- 
house Electric and Manufacturing Company. 
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At the conclusion of the afternoon session there was an inspection of the 
Physics Laboratories of Hunter College. The inspection was followed by the 
dinner of the Metropolitan Section in the restaurant of Hunter College. 


Evening Session 
Invited Papers: 


The New Home of the American Institute of Physics. H. A. BARTON, American Institute of 


Physics. 


The Chemical Age of Plenty. F. C. WHiTMorE, Pennsylvania State College. 


Absorption of -Radiation in Terrestrial Materials. 
C. W. SHEPPARD, Massachusetts Institute of Technology.— 
As part of a research program on the role of radioactivity 
in petroleum genesis calculations have been made of the 
absorption of y-rays in terrestrial materials. The intensity 
of ionization has been calculated for a steel ionization 
chamber having thick walls and containing air at one 
atmosphere. In the case where the chamber is situated 
above a semi-infinite solid absorber containing mpra, mTh, 
and mx grams of Ra, Th, and K per gram, respectively, 
the number of ion pairs per cc per second will be 
1.4X 10"mpat+2.4X10'm7,+27mxK. The expression was 
tested by inserting the measured Ra content! and Th 


W. S. GorTON 
Secretary-Treasurer 


content? of selected samples from the Millikan rock series, 
together with an estimate of the K content. The results 
were compared with Millikan’s measured values using the 
factor (1/13.8) to reduce the readings to 1 atmosphere.' 
No systematic disagreement was found. Large individual 
deviations in both directions indicated sampling irregu- 
larities. The lower value of ionization estimated by Hess* 
may be attributed to under-estimation of the effect of 
scattered low energy quanta on the ionization chamber. 

1R. D. Evans and R. W. Raitt, Phys. Rev. 48, 171 (1935). 

2R. W. Raitt, Ph.D. Thesis, California Institute of Technology 


(1935). 
3V. F. Hess, Terr. Mag. 46, 409 (1941). 





Proceedings of the American Physical Society 


MINUTES OF THE MEETING OF THE NEW ENGLAND SECTION HELD AT THE CONNECTICUT COLLEGE 
FOR WoMEN, NEw LONDON, CONNECTICUT, OCTOBER 16, 1943 


‘HE twenty-second regular meeting of the 
New England Section of the American 
Physical Society was held at the Connecticut 
College for Women, New London, Connecticut, 
on Saturday, October 16, 1943. Forty-five mem- 
bers were in attendance. The invited papers were 
as follows: 


Discussion on “Problems Connected with the V-12 Pro- 
gram” led by V. E. Eaton, Wesleyan University. 

Physics at Connecticut College. G. K. DAGHLIAN, Con- 
necticut College. 

Some Crude, Simple, Large Scale Models of Apparatus 
Developed for Teaching First-Year College Physics. 
J. Barton Hoac, U.S. Coast Guard Academy. 

Fast Neutron Absorption in Gases, Walls, and Tissue. 
Griapys A. ANsLow, Smith College. 

Is Hardness a Physical Property of Solid Matter? S. R. 
WitiiaMs, Amherst College. 


But one ten-minute paper was presented for 
which the abstract follows: 


Application of Dirac’s New Field Theory to the Calcula- 
tion of Higher Order Matrix Elements. Sauicui Kusaka, 
Smith College.—Dirac’s new method of eliminating the 
divergences in field theories is applied to the calculation 
of higher order matrix elements which diverge in the usual 


theory. For an electromagnetic field interacting with 
charged particles, all these matrix elements vanish, and 
hence this theory coincides with the theory of radiation 
damping developed by Heitler, Wilson, and others in 
which the matrix elements which diverge in the usual 
theory are neglected. For all types of meson fields inter- 
acting with nucleons, however, Dirac’s theory gives finite 
but non-vanishing values for the higher order matrix 
elements, and hence does not agree with the current theory 
of radiation damping. Moreover, these matrix elements are 
of the same order of magnitude as the lowest order matrix 
element even for non-relativistic energies, and become 
much larger for higher energies. Thus the perturbation 
theory is inadequate for the treatment of these problems. 


At the business meeting a resolution was pre- 
sented to the meeting in commemoration of 
Arthur P. R. Wadlund of Trinity College, who 
was Vice Chairman of the Section at the time of 
his lamented death. The following officers were 
elected for the ensuing year: 

Chairman: Gladys A. Anslow, 

Vice Chairman: Morton Masius, 

Secretary-Treasurer: Mildred Allen, 

Members of Programme Committee: Clarence 
E. Bennett, W. W. Stifler. 

MILDRED ALLEN, Secretary-Treasurer 
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